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OSCILLATING TEMPERATURE CAPILLARY ELECTROPHORESIS 

AND USES THEREFOR 

RELATED APPLICATION 

This application claims the benefit of the U.S. Provisional Application No. 

5 60/ , entitled, "Oscillating Temperature Capillary Electrophoresis and Uses 

Therefor," by Per Olaf Ekstrom, filed November 12, 2003, having attomey docket 
number 3637.1000-000. 

The entire teachings of the above application are incorporated herein by 
reference. 

1 0 BACKGROUND OF THE INVENTION 

In 1953, Watson and Crick hypothesized the stmctm-e of the DNA molecule. 
Approximately five decades later the draft sequence of the human genome was 
completed. Decoding of the human DNA sequence through the Hmnan Genome 
(HUGO) project differences in the genomes were analyzed. The most firequent DNA 

15 variations, have elucidated single nucleotide polymorphisms (SNPs), are considered 
to be useful polymorphic markers for genetic studies of pharmacogenetics and 
polygenic traits. A worldwide effort has achieved an accumulation of millions of 
SNPs stored in public databases. However, most of these SNPs have been identified 
by examination of a limited mmiber of individuals, and infomiation on their allele 

20 firequencies is lacking or tentative. Furthemaore, studies have shown that allele 
fi-equencies vary widely between different ethnic populations. Thus validation of 
SNPs and estimation of their allele firequencies, especially for each ethnic group, are 
required before these markers can be used for genetic studies. 

Current efforts in discovery and screening of DNA variations ^oint 

25 mutations, SNPs and other polymorpliisms) are expected to yield valuable 

infonnation on potential genetic risk factors as well as information on genetic 
variants that could lead to enhanced susceptibiUt}^ for certain diseases. It is expected 
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that discovery and massive screening of DNA polymorphisms will become essential 
for taylor-made drugs as well as disease gene association studies. Screening for and 
identification of SNPs and other polymorphisms will therefore require methods for. 
separating variant DNA sequences fi'om large samples. 

5 SUMMARY OF THE INVENTION 

Described herein is a new type of DNA variation (polymorphism, mutation 
or SNP) screening technology based on capillary array electrophoresis. Sample 
DNA jfragments of known sequence are PGR ampUfied and detected based on their 
differential migration in a polymer-filled capillary. A cycling (oscillating) 

10 temperature is used to compensate for local fluctuations in temperatures across the 
multi-capillary array. In addition, the application of short periodic temperature 
cycles allows a user to continue injecting subsequent sample plates between the 
temperature gradient cycles before the earlier samples appear in the detector. Using 
this novel approach, a dramatic and unexpected increase in separation throughput by 

15 more effective utihzation of the separation capacity (volume) of each capillary is 
demonstrated. 

In one embodiment, the present invention is directed to a metliod for 
separating nucleic acids comprising electrophoresing a sample applied to a gel 
electrophoresis matrix in a capillary, wherein during electrophoresis, the temperature 

20 of the matrix is cycled at least two times between a high and low temperature. In a 
particular embodiment, the nucleic acids to be separated are DNA firagments 
comprising one or more polymorphic sites, hi another embodiment, allelic variants 
at the one or more polymorphic sites are separated. In a particular embodiment, the 
temperature is initially at a high temperature and the first cycle is fi*om a high 

25 temperature to a low temperatm*e. In another embodiment, the high temperature 
and/or low temperature is different during successive cycles, hi another 
embodiment, the temperature is cycled fi-om about 2 to 60 times. In a particular 
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embodiment, tlie temperature is cycled about 20 times. In another embodiment, the 
high temperature is about 3 °C higher than tlie low temperature. In another 
embodiment, the temperature is between about 2 °C and about 15 °C higher than the 
lower temperature. In another embodiment, the higher temperature is between about 
5 3 °C and about 10 °C higher than the lower temperature. In one embodiment, the 
high temperature is less than about 80 In a particular embodiment, the low 
temperature is about 40 In one embodiment, the high temperature is between 50 
°C and 75 °C. In another embdoiment, the low temperature is between 40 °C and 50 
°C. Ina particular embodiment, the method further comprises detecting dsDNA after 

10 electrophoresis. In a particular embodiment, after the desired number of temperature 
cycles have been completed, the temperatui-e of the gel matrix is such that DNA 
remains double-stranded. In another embodiment, the temperature oscillations are 
ramped to provide optimal separation of the alleles. 

In another embodiment, the present invention, is directed to a method for 

1 5 estimating allele frequency comprisirig: electrophoresing a sample applied to a 

capillary gel electrophoresis matrix, wherein during electrophoresis, the temperature 
of the matrix is cycled at least two times, wherein one cycle is from a high 
temperature to a low temperature or fi'om a low temperature to a high temperature, 
thereby separating DNA molecules in the sample; and quantifying the variant 

20 sequences of the separated DNA molecules thereby providing an estimate of the 
allele frequency for each variant DNA molecule. In a particular embodiment, the 
method further comprises detecting dsDNA after electrophoresis. In one 
embodiment, a&&r the desired nmnber of temperature cycles have been completed, 
the temperature of the gel matrix is such that DNA remains double-stranded. 

25 In yet another embodiment, the invention is directed to a method for 

detecting a microhaplotype comprising separating DNA fragments comprising a 
sequence comprising two or more polymorphic sites of the microhaplotype, wherein 
the fragments are separated by capillary electrophoresis performed with two or more 
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temperature oscillations between a high and a low temperature. In a particular 
embodiment, the temperature oscillations are ramped to provide optimal separation 
of the microhaplotype. 

BRIEF DESCRIPTION OF THE DRAWINGS 
5 FIGS. lA-lC are graphs showing temperature over time during different CE 

methods. FIG. 1 A shows multiple-injection CE. FIG. IB shows multiple- injection 
TGCE. FIG. IC shows multiple-injection CGCE. The upper drawing illustrates the 
temperature profile. The simulated positions of separated sample zones are at the 
bottom. 

10 FIG. 2 is a spectral profile showing a comparison of a two temperature 

profile mode. Single samples analyzed for single nucleotide polymorphisms in the 
LTA gene by temperature gradient capillary electrophoresis. Genotypes of the 
sample are heterozygous GA; thus four peaks are identified as GC, AT, AC and GT, 
respectively. Lower part of the figure display progranmied temperature profile of 

15 the analysis chamber during electrophoresis. 

FIGS. 3 A and 3B show experimental measurements of temperature profiles 
inside MegaBACE™ 1000 capillary chamber. FIG. 3 A depicts the typical profile of 
a single-sweep gradient mode. FIG. 3B shows a cycling temperature gradient mode. 
The MegaBACE™ lOOO used here was a modified "HOT" MegaBACE™ at 

20 Molecular Dynamics. The standard temperature limit, 50 *=*C, on the MegaBACE™ 
was raised to 65 °C. 

FIG. 4 is a fluorescence and temperature profile depicting the analysis of 
several samples for APC mutation detection using multiple-injections and fast 
cycling temperature gradient conditions. The three injections were performed every 

25 fourth temperature cycle. First injection is C homozygote, second injection is T 
homozygote and the last injection is C/T heterozygote. 
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FIG. 5 is an illustration of liigh throughput SNP scoring by CTCE in multi- 
capillary format. Five different sample plates (FAM-labeled PGR fragments of 
BRCA2) were consecutively injected during a cycling temperature gradient. Allele 
scoring was performed against a TMR-labeled internal mutant showing heterozygote 
5 profile with both allele peaks. The intemal standard was also used as an "injection 
marker" for consecutive data processing by SNP Profiler software. 

FIG. 6 is a schematic diagram depicting an overview of CTCE sample 
processing workflow. 

FIG. 7 is a chart showing the relative thermodynamic stability of mutations 
10 as compared to wild-type fragments. Wild-type KRAS exon 1 with the sequence 
GGT and GGC in codon 12 and 13, respectively, is coded as 100% and the scale 
indicates relative difference in tlie thermodynamics profile of the mutants. All G to 
C mutations proved to be more stable than wild-type, and will thus separate out 
before wild-type peaks in electropherograms, except for a GGC to GCC mutation 
15 that was only slightly more stable than wild-type, and will comigrate with wild-type. 
All other mutations were less stable than wild-type and will separate after wild-type 
peaks. Note that some of the mutations had veiy similar thermodynamics. 

FIG. 8 is a representation of electropherograms of all oncogenic mutations in 
KRAS exon 1 and standard analyzed with CTCE. Separation between the wild-type 
20 homoduplex and the heterozygous peaks was achieved in all mutated samples. 

Separation between the wild-type homoduplex and mutant homoduplex was foimd 
in 1 1 out of 12 mutants. The mutant homozygous peak of GGC to GCC 
transversion in codon 13 comigrated with wild-type at the selected denaturing 
conditions. This observation corresponded with the theoretical melting properties of 
25 the GCC mutation as compared to wild-tjpe sequence. The thin line represents the 
sample and the bold line the standard. The peaks of the standard are named S1-S5 
from left to right (SI, mutant homoduplex from CGT; S2, mutant homoduplex from 
TGT; S3, mutant homoduplex from AGC; S4, mutant heteroduplexes from TGT and 
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CGT; S5, mutant heteroduplex jfrom AGC). The peaks of the sample are named WT 
(wild-type), M (mutant), HI (heteroduplex with shortest migration time) and H2 
(heteroduplex with longest migration time). Note that some of the samples contain a 
single-strand peak produced 
5 vmder the PGR amplification or under the CTCE analysis. 

FIG. 9 is a diagram showing the Genetic Profiler analyzed migration times of 
all peaks in the electropherograms, and the relative migration times between peaks 
were calculated by Excel. Individual mutations proved to have distinct migration 
times when compared to the peaks of the standard and other mutations. Migration 

10 times that were able to separate different mutations were selected and a classification 
tree was made. For example, a TGT mutation was characterized first by a shorter 
migration time of HI than S4, and migration time of S4 minus migration time of HI 
was less than 12. Finally, the migration time of the mutant homoduplex (M) was 
longer than for the wild-type peak (WT). No other mutations had these 

15 characteristics. 

FIG. 10 is a graph showing electropherograms of 5 samples heterozygous in 
TNF target single-nucleotide polymorphism (SNP). Ratio between the G and A 
alleles was calculated by measuring tlie area under the peaks and foimd to be 0.996 
(one SD = 0.028). 

20 FIG. 1 1 is a graph showing electropherograms of the 4 pools representing 

2000 alleles each, analyzed for the single-nucleotide polymorphism (SNP) in TNF 
by denaturant capillary electrophoresis (DCE). The area under the peaks was 
measured and the allele frequencies calculated. Pool 3 revealed a large single-strand 
pealc (marked X) created either during the PGR or during electrophoresis. The 
25 single strand did not mterfere with the measurement of pealc areas. 

FIG. 12 is a table showing allele and genotype frequencies of anal)/zed SNPs. 
FIG. 13 is a plot showing electropherograms of three samples with different 
genotypes. The samples were analyzed by denaturant capillary electrophoresis 
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(DCE) with internal Tamra-labeled standard (bold line). Upper part shows the 

n 

I. 

homozygous GG genotype, the middle part shows tlie genotype AA, and lower part 
shows the heterozygous GA genotype. 

FIG. 14 is a plot showing five heterozygous samples injected in 8.5 minute 
5 intervals during denaturant capillary electrophoresis (DCE) with a gradient of 54 °C 
to 50 °C in 1 °C steps. Each temperature was held for 1 minute. The gradient 
was repeated 8 times. 

FIG. 15 A-E are electropherograms from cycling temperature capillary 
electrophoresis of five genotypes representing the five microhaplotypes combination 
10 found. Sample alleles are depicted with thin line and the intemal standard are 

presented as bold line. Microhaplotypes are given as the three polymorphic bases in 
each allele. F or example, the first peak in the intemal standard is IVS38-8 T 5557 
G 5558 A, is labeled as TGA. Note that heteroduplexes are not shown in the 
electropherograms. 

1 5 FIG. 16 is a table showing the total nimibers of haplotypes observed in the 

normal and rectal cancer groups. Normalized melting behavior of microhaplotypes 
are presented under calculations. High temperature property is indicated by negative 
values for the fragments. For positive values, the fragment would be less stable as 
compared to the wild-type and consequently elute after the wild-type. 

20 FIG. 17 is a table showing the observed genotypes in the normal population 

(n = 3,526) and rectal cancer patients (n = 151) and allele frequencies of each SNP. 

FIG. 18 is a graph showing relative differences (y-axis) in stabiUty between 
G and A alleles in lymphotoxin alpha at different melting probabilities (x-axis) are 
shown. The thermodynamics calculated by the use of the Poland web page (H) 

25 predicted the observed stability difference betv^^een the alleles (9) at 50% melting 
probability. WioMelt (F) predicted difference in separation at all probabihty levels. 
At 50% melting probability, a good theoretical prediction of observed 
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thermodynamic differences was obtained with both the Poland home page and 
WinMelt. This level of probabiUty was used for fixrfher analyses, 

FIG. 19 is a plot showing the correlation between predicted thermodynamics 
at 50% melting probability with WinMeh and relative migration time in DCE under 
5 optimal denaturing conditions. A linear correlation with r^ = 0.97 was found. 

FIG. 20 is a representation showing temperature measured at three different 
locations with calibrated external temperature sensors in the capillary chamber at 
three different temperature settings. Temperatures in different parts of the capillary 
chamber were recorded continuously for a period of 15 min. The dotted line 

10 represents the intended temperature in the chamber. Temperaturesmeasured around 
the capillaries in the front are indicated with an A, the capillaries in the middle with 
a B, and the capillaries at the back with a C. The upper part of the fiigure shows the 
temperature recorded at different locations with a preset constant temperature in the 
chamber. In the middle part of the figure a gradient going through the theoretical 

15 denaturing temperature is employed. Li the lower part of the figure, cycling of 
temperature around the intended temperature is illustrated. 

FIG. 21 is a plot showing relative retention of the two homoduplexes and 
two heteroduplexes in a heterozygous sample analyzed with DCE at different 
predefined denaturing temperature. One standard deviation is noted in the figure for 

20 each homoduplex (n = 1 1) and heteroduplex (n = 1 1). 

FIGS. 22A and 22B are plots showing electropherograms. In FIG. 22A, two 
electropherograms of the same PGR firagment, analyzed in the same run and in 
differentcapillaries under constant denaturant conditions, are shown. A difference in 
separation is observed between the two capillaries. In FIG. 22B the same PGR 

25 fragments are analyzed with cycling temperature, witli a gradient from 48 °C to 52 
°C and cycled 15 times. Each temperature was held for 30 s. With these denaturing 
conditions, equal migration times were observed in different capillaries. 
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FIG. 23 is a plot showing a KRAS exon 1 internal standai-d with three 
homoduplex peaks and two heteroduplex peaks was analyzed with CTCE in 96 
capillary format. The peaks were measured in 649 individual samples and presented 
as box plots to demonstrate the reproducibihty in DCE analysis with cychng 
5 temperature. 

FIG. 24 is a theoretical melting profile of twelve double-stranded DNA 
fiagments as analyzed by the WinMelt computer program. The following DNA 
sequences are depicted: 1) BRCA 1, 2) BRCA 1, 3) MTHFR, 4) OPSIN, 5) 
MTHFR, 6) MTHFR, 7) CBS, 8) NQOl, 9) DPYD, 10) DPYD, 1 1) DPYD, 12) 
10 CTLA-4. Theoretical meltuig temperatm*e is depicted on the y-axis, and number of 
base pairs on the x-axis. Note that the maximmn difference in thermodynamics 
between two fragments is 9.4 K. 

FIG. 25 is a plot showing electropherograms of all twelve fragments when 
analyzed with CTCE and cycling denaturant temperature from 59-47 ^'C and 5 
15 cycles. The fragments are arranged as in FIG. 24. Note the baseline separation 
between the two homoduplexes and between homoduplexes and heteroduplexes in 
all electropherograms. 

FIG. 26 is a schematic showing the difference in migration times of all peaks 
after CTCE analysis with cycling denaturant temperature from 59-47 °C and 5 
20 cycles. The results are illustrated as average (n = S) peak maximum with one 
standard error of the mean. 

FIG. 27 is a plot showing temperature measured at three different locations 
with calibrated extemal temperature sensors recorded continuously during 
electrophoresis. The cycling condition was temperatures programmed from 59-47 
25 ''C for 5 cycles. Temperature measured in the front, middle and back of the chamber 
are as indicated in the legend. 

FIG. 28 is a graph illustrating the difference in migration time between the 
peaks of a heterozygous PCR sample at different denaturing conditions. Each point 
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in the graph represent separation at 3 K gradient with average gradient temperature 
starting at 42.5 °C and increasing by one °C. Note that separation of individual 
peaks can be achieved within a quite large temperature frame. 

FIG. 29 is a plot showing observed denaturing temperature in the capillary 
5 chamber at three different locations. The programmed denaturing temperature was 
from 47 **C to 55 and with 7.5 cycles. Temperature measured in the front, middle 
and back of the chamber are as indicated in the legend. 

FIG. 30 is a plot showing the difference in migration times of all peaks after 
CTCE analysis with cycling denaturant temperature from 47-55 "^C and 7.5 cycles. 
10 The results illustrate average peak maximum and one standard error of tlie mean. 
Note that fragments with high theoretical thennodynamics (number 6 and 7) have 
insufficient separation. 

FIG. 31 is a plot showing the theoretical thermodynamics of a target 
sequence. Base substitution A->G in the lower domain changes the melting 
15 properties as shown in the enlargement 

FIG. 32 is a plot showing three samples analyzed for mutation in KRAS 
exon 1 . The upper electroplireogram shows the wild-type sample, whereas the two 
lower samples harbor GGT>GTT mutations. Peaks numbered 1 to 4 represent the 
wild-type homoduplex, mutant homoduplex and heteroduplexes, respectively. Note 
20 the different mutant fraction in the two lower electropherograms. 

FIG. 33 is a plot showing analysis of SNP in IL-4 gene by DCE. 
Identification of genotypes by co-elution with an internal standard (bold line) made 
up of both alleles. The genotypes from top to bottom are TT, CT and CC, 
respectively. Heteroduplexes are not shown. 
25 FIG. 34 is a table showing microhaplotype combination for two SNPs 

separated by 57 bp. 
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FIG. 35 is a plot showing five samples inicrohaplotyped according to an 
internal standard (bold line). Genotypes are given as allele combination as 
numbered in FIG. 34. Heteroduplexes are not shown. 

FIG. 36 is a plot showing quantitative analysis of gene coijy nxraibers in 
5 Opsin gene array. The ratio of L/M based on area under the peaks is 1/3. Note the 
non-template, polymerase-mediated adenine addition (a), which does not interfere 
with the measurements. 

FIG. 37 is a plot showing allelic imbalance analyzed by DCE. The ratio of 
the two alleles in the upper sample is 1 : 1 as expected for a hetero2ygous sample. In 
10 tlie lower electropherogram a clear imbalance between alleles can be seen. 

DETAILED DESCRIPTION OF THE INVENTION 

Described herein is an oscillating temperatme gradient electrophoresis 
system useful for separating biomolecules, e,g.^ DNA. The methods described 
herein can be used in combination with multiple-injection technology to 

15 dramatically increase the sample throughput (Mansfield, E. et al, 2000, U.S. Patent 
No. 6,156,178). The present work extends technology for detection of DNA variants 
into a high tliroughput mode. Methods for screening multiple samples for the 
presence of somatic mutations after transferring a liigh-resolution technique of 
CDCE fi-om single-capillaiy fomiat to a commercial capillary array instrument had 

20 been described (Bjorheim, J. et al, 2002, Anal Biochem,, 304:200-205). Additional 
advantages of applying an oscillating temporal temperature gradient rather than 
maintaining an accurate constant temperature for increased reproducibility in a 
conmiercial multi-capillary instrument were later described (Minarik, M. et al, 
2001, poster presented at ICHG Conference, Vieima, Austria). 

25 Samples suitable for the present invention can be biological samples 

obtained firom, for example, animals, plants, viruses or bacteria. Animal samples 
can be, for example, manmialian in origin and derived from, for example, blood. 
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seruin, cells tissue. For example, a biological sample from a test subject (a "test 
sample") containing a biomolecule of interest, e.g., genomic DNA, RNA, protein, 
polypeptide or cDNA, is obtained from an individual or several samples can be 
collected from multiple individuals in a population. The individual can be an adxilt, 
5 child, or fetus. The test sample can be from any source that contains the desired 
biomolecule, preferably DNA, such as a blood sample, sample of amniotic fluid, 
sample of cerebrospinal fluid, or tissue sample from skin, muscle, buccal or 
conjunctival mucosa, placenta, gastrointestinal tract or other organs. A test sample 
of DNA from fetal cells or tissue can be obtained by appropriate methods, such as by 

10 anmiocentesis or chorionic villus sampling. 

Among the various experimental techniques for detection of DNA variations, 
DNA sequencing is the most imiversal approach. Although sequencing has been the 
"gold standard" for detecting polymorphisms for many years, a variety of alternative 
approaches for detection and screening have emerged, overcoming the high cost and 

1 5 occasional lack of sensitivity. These methods can rely on hybridization, 
allele-specific enzymatic reactions such as PGR, minisequencing, strand 
displacement, and/or cleavage, or on differences in physicochemical properties of 
the DNA variants such as the melting equilibrium, or affinity towards stationary 
phase. The results of these techniques can be read directly from the reaction mixture 

20 using plate readers or following separation of nucleic acids by capillary 

electrophoresis (CE), HPLC or mass spectrometry. The techniques employing CE, 
HPLC or mass spectrometry offer an additional means of identification of the 
variation from a characteristic pattem of migration, elution or mass spectra. With 
many of the above-mentioned screening techniques, the sample throughput can be 

25 increased by processing multiple samples in parallel. This feature becomes useful 
when confronted with the vast amount of potential sequence targets to be 
scmtinized. 
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Electrophoresis has traditionally been used for separation of DNA fragments, 
DNA sequencing as well as general fragment analysis. A large family of slab-gel 
techniques dedicated to detection of DNA variations is based on differential melting 
of wild-type (taken to be a particular reference sequence) and mutant (a sequence 
5 that varies from the wild-type sequence) DNA fragments translated into an 

observable retention difference through electrophoretic sievuig. One technique that 
utilizes this principle is denaturing gradient gel electrophoresis (DGGE), in which 
amplified fragments of wild-type and mutant sequences are resolved during their 
migration in a slab gel containing a gradient of chemical denaturant. DGGE is well 

10 estabhshed in clinical diagnostics due to its relative simplicity and ability to resolve 
close to 100% of mutations present in a sample for a given target sequence. 
Following DGGE, other variants of slab-gel mutant separation methods were 
developed including temporal temperature gradient gel electrophoresis (TTGE), in 
which the temperature is varied throughout separation and constant denaturing gel 

1 5 electrophoresis (CDGE), where the separation takes place at predetermined constant 
denaturing conditions. 

In addition to the electrophoretic methods, an approach based on denaturing 
HPLC (dHPLC) has been used to identify polymorphsims. dHPLC uses an ion-pair 
chromatography separation principle combined with precise control of the colxmm 

20 temperature and optimized mobile phase gradient for separation of mutant 

heteroduplexes. dHPLC can be easily automated and offers an option to collect the 
isolated heteroduplexes for further identification or confirmation by sequencing. 
However, the main potential of DGGE and dHPLC is mostly in discovery of novel 
mutations rather than screening due to their relatively low throughput of 

25 approximately 5 nun per analysis. 

The transition in DNA separations from traditional slab-gel electrophoresis 
to CE systems started in early 1990 and was later ftirther accelerated by the Human 
Genome Project. Separation of heteroduplexes was achieved at different 



wo 2005/047541 



PCT/IB2003/005304 



-14- 

temperatures controlled by means of Joule heating through adjustment of the 
separation voltage, hi 1994, a metliod referred to as constant denaturing capillary 
electrophoresis (CDCE) was introduced as a capillary analog to CDGE, a constant 
denaturing electrophoresis performed on slab gels. For CDCE, the separation is 
5 carried out at an accurately maintained constant temperature at which the homo- and 
heteroduplex forms of wild-type and mutant sequences exhibit the best separation. 
An alternative approach was later used that overcomes the requirement for a very 
accurate temperature control by applying a simple temporal temperature gradient. 
Applying a temperature gradient, rather than a constant temperature, is useful 
1 0 especially in multi-capillary systems where maintaining accurate temperature across 
all capillaries is difficult With a temperature gradient, each capillary reaches its 
temperature optimum, even if there is a difference in absolute temperature values 
among capillaries. The duration of such optimum separation conditions only 
depends on the overall slope of the temperature gradient. 
1 5 The temperature gradient approach can be appHed to many existing 

commercial multi-capillary CE systems with no additional requirements on 
instrument hardware. The typical run time is less than 60 minutes, depending on the 
resolution requirement (i,e,, gradient slope). 

The present invention encompasses methods to improve separation of, for 
20 example, different alleles of a DNA fragment. An "allele" is a particular sequence 
version present at a "polymorphic site", i.e., a site that, in a particular population of 
molecules, can have more than one nucleotide present. For example, a specific 
position of a DNA fragment of a human chromosome might have an adenine in one 
molecule and a cytidine in another molecule at the same position. The position 
25 would be a "polymorphic site" having at least two possible alleles: the "A" allele 
and the "C" allele. Since the fragments might otherwise be identical in sequence 
(see Example 7 for a discussion of "microhaplotypes"), the problem of separating 
the A allele from the C allele is significant. Described herein are CE methods for 
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separating DNA fragments that might differ in sequence at only one position. 
Samples are PCR-amplified to attach a "GC-clamp" (a sequence high in guanine and 
cytidine residues), and subsequently separated by Capillary electrophoresis. By 
allowing the temperature to oscillate during the electrophoresis run, it was 
5 unexpectedly observed that separation of biomolecules was dramatically improved. 

The methods of the present invention can be used in conjunction with known 
and yet to be discovered electrophoresis methods to improve, for example, 
throughput. For example, temperature oscillations can be introduced dimng 
multiple injection capillary electrophoresis to further improve resolution and 

10 increase the amount of sample that can be separated during a given time period. 
Additionally, the invention can be used in parallel in multi-capillary systems. 

In a multiple-injection experiment, samples are serially injected into the 
capillary (or array of capillaries) in periodical time intervals sepai-ated by short 
application of separation voltage. This approach is an efficient usage of the 

15 separation capacity (migration volume) of the capillary colimm. Following injection 
of the first sample (or a set of samples in capillary array), a so-called interval voltage 
is applied for a sufficient period of time (typically 2-5 min for short 
oUgonucleotides) preventing overlap of the slowest peak from the sample with the 
fastest migrating peak from the next sample. After this period, the next sample is 

20 injected resulting in a continuing process during which the first set of samples is 
reaching detectors, the following set is separating and new ones are being injected. 
In the most common version, Ibe injections are repeated several times before the 
fastest migrating peaks from the first sample reach the detector. After the final 
injection, a run voltage is applied for a longer period of time (30-60 minutes) to 

25 drive all peaks from all injections pass the detection window. The maximimi 
number of injections are determined experimentally. 

Most commercial capillary array instruments allow controlling temperature 
during the run. FIG. 3 shows the temperature profile of a single-sweep gradient (A) 
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and atypical cycling gradient (B) recorded inside the capillary chamber of 
MegaBACE 1000 capillary array instrument equipped with high temperature setting. 
The unchanged separation performance of the cycling gradient compared to a 
single-sweep gradient allows for the application of the multiple-injection method. 
5 Depending on the actual application, the samples can be periodically injected either 
in every gradient cycle or once in every few cycles (e.g., every cycle, every 2 cycles, 
every 3 cycles, every 4 cycles, every 5 to 10 cycles, every 7 to 20 cycles, every 15 to 
40 cycles, etc.). A single cycle can start at either high temperature limit or the low 
temperature limit. An example of multiple-injection cycling temperature gradient 
. 10 capillary electrophoresis (CGCE) with fast temperature cycling (5 minutes per cycle) 
is shown in FIG. 4. APC mutant samples were injected in every fourth gradient 
cycle. During each cycle, a temperature gradient going from 55 °C to 53 ""C and 
back to 55 ^'C was applied. An intemal standard containing 1 : 1 mixture of 
individually amplified wild-type and mutant fragments labeled with a TMR dye was 
15 included in each sample well. Since no heteroduplexes were formed following the 
PGR, the resulting separation pattems consist only of homoduplexes. It can be seen 
that a fall separation of both wild-type and mutant homoduplexes of the intemal 
standard was achieved within the individual interval windows. Clearly, a very high 
resolution of the peak separation is achieved imder the cycling temperature gradient 
20 conditions, allowing identification of both alleles in a heterozygous sample by 

comparison to the peak pattem of the standard. Apart from just detecting a presence 
of DNA variation from the characteristic peak pattem in heterozygous samples, the 
complete separation of the two homoduplexes allows direct identification of 
homozygous genotypes. This is the key in automated SNP genotyping, where the 
25 two homozygous genotypes can be directly scored. An example of high-throughput 
analysis and automated scoring of SlsfPs is shown in FIG. 5. A BRCA2 SNP was 
screened in various patients using the five injections with 45 °C - 43 °C cycling 
temperature gradient. Each sample included a 6 - carboxy - N,N,N%N' - 
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tetramethylrhodamine (TAMRA)-labeled internal standard. The data was processed 
using SNP Profiler software which allows assignment of individual injection 
windows. The genotypes are directly determined from the positions of mutant 
homoduplex peaks co-eluting with the internal standard homoduplex peaks. 
5 The total runtime of the experiment shown in FIG. 5 was less than 2 hours. 

It was estimated that, on a common 96-capillary instrument, five injections could 
feasibly be performed without any adverse effect on the separation matrix. The total 
runtime would then be 20 min of initial "dead" volume + 5x15 minutes separation 
window +15 minutes final electrophoresis =110 min. Considering 10 minute 

10 additional periods required for gel replacement in betv^^een runs, a total of 12 of 
similar multi-injection runs can be performed in less than 24 hoxxrs of operation. 
This represents an overall throu^put of 5760 samples in 24 hours. Commercial 
genetic analyzers are usually equipped with four dye channels. Considering that up 
to three fluorescent dye chaimels can be used to detect unknown samples (the last 

1 5 channel is assigned for the internal standard), it is clear that the potential capacity 
can be fiirther increased 3-fold up to 17,200 samples in 24 hours on a single 
96-capillary instrument. Unlike in otlier mutant or SNP scoring methods, the 
presented technolog)' includes a very straightforward workflow shown in FIG. 6. 
Following the original PGR amplification, there is no sample cleanup required. 

20 Using a 96-format, PCR themiocycler enables complete automation firom sample 
preparation to the multiple-injection CE analysis using robotic plate handlers. 

The methods of the present invention can be used in conjunction for 
separating biomolecules (DNA, RNA, DNA:RNA hybrids, polypeptides, proteins, 
carbohydrates, etc.) that are conventionally separated by electrophoresis. The 

25 methods described herein apply to those biomolecules that are affected by 

temperature variation during separation. The method is not limited to CE, but can 
be applied to other electrophoresis methods as well. For example, temperature 
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oscillations can be introduced during slab gel electrophoresis, polyacrylamide gel 
electrophoresis (PAGE), sodium dodecyl sulfate PAGE, etc. 

The temperature oscillations can be performed with a high temperature 
boundary approximating the temperature at which the GC-clamp melts. For 
5 example, a GC-clamp having a melting temperature of 85 °C will have a practical 
high temperature limit of 85 °C during electrophoresis. The practical lower 
temperature limit for temperature oscillations during electrophoresis is the minimal 
melting temperature of AT-rich DNA sequences. In practical terms, the lower 
temperature limit is about 35 °C. Thus, temperature oscillations can be perfonned 

10 from about 35 °C to about 85 °C, depending on the fragment to be separated and the 
resolution desired. The difference between the high and the low temperature Umits 
for the oscillations can be, for example, about 1 "^C, about 2 "^C, about 3 °C, about 5 
°C, about 7 °C, about 10 °C, about 15 °C, about 20 °C, about 25 °C, about 30 °C, 
about 35 °C, about 40 °C, about 45 ^'C or about 50 °C. Additionally, the high and/or 

1 5 low temperature limits for the oscillations can vary from cycle to cycle. For 

example, the first cycle can have a high temperature T^i and a low temperature Tlj, 
and have a different high temperature during the second cycle, Th2 and T^. Each 
cycle can have a different high and low temperature, or they can be the same for the 
purposes of the present invention. 

20 One of skill in the art would know how to vary electrophoresis conditions, 

e.g., through the addition of compounds such as, for example, urea, such that DNA 
melting temperature will be affected. Conditions can be estabUshed to generally 
lower the melting temperature of dsDNA, for example, by adding urea to the gel 
and/or electrophoresis buffer, or conditions can be altered to lower or raise melting 

25 temperatures of specific dsDNA sequences. 

Additionally, the rate at which the temperature cycles from the high 
temperature to the low temperature ("ramping" of temperature), can be altered to 
affect separation. For example, if two fragments to be separated differ in melting 
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temperature by 1 °C, a longer ramp time would improve separation. Separation will 
depend on the amount of time one fragment is partially denatures while the other 
fragment is double-stranded. By extending the ramp time of oscillations, the time 
during electrophoresis spent between the two melting temperatures is increased, 
5 thereby improving separation. 

As a result of the invention described herein, highly sensitive and accurate 
methods of detecting polymorphic DNA and detemiining allele frequencies are now 
available. 

EXEMPLIFICATION 

10 Example L Cycling gi'adient capillaiy electrophoresis: A loM^-cost tool for 
high-throughput analysis of genetic variations. 

In the present work, the general principle that applying a temporal 
temperature gradient in CE can further be extended into applying a periodical 
temperature cycle is demonstrated. PGR is employed with one of the primers 

1 5 extended by a high-melting domain (the "GC-clamp") to amplify a target DNA 

sequence surrounding the mutant or SNP marker. The PGR conditions are specific 
for each tai'get sequence. 

The application of periodic cycles allows a better compensation of the local 
temperature fluctuations inside the multi-capillary oven. Rapid gradient cycling 

20 with rates of up to several cycles per minute showed better results in comparison 
with slower cycling intervals. In addition, the instrument hardware does not appear 
to be able to follow the rapid temperature changes inside the chamber resulting in a 
relative constant average temperature of the optical components. Periodic 
appUcation of tlie temperature gradient enables usage of multiple injection 

25 technology, in wliich different samples are serially injected between the cycles and 
separated under the same revolving temperature conditions. Multiple-injections 
allow for significant increase in sample throughput. Gonsidering the ease of sample 
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preparation, PGR directly followed electrophoresis without any post-PCR treatment 
such as desalting or removal of unincorporated primers. Evaluation of the mutant 
presence or SNP genotype is performed solely based on an intemal standard running 
in a separate spectral channel. In situations where a slower migrating PGR fragment 
5 would directly overlap with a faster migrating primer peak from a subsequent 

injection, the spacing between the injections was adjusted. For a given mutant/SNP 
marker, the peak distance is very reproducible, therefore the injection spacing can be 
optimized accordingly. This technique represents a cost effective, simple and 
powerful tool for high- throughput scoring of DNA mutants and SNPs. 

10 Ghemicals 

All experiments were performed using standard MegaBAGE. buffers and 
MegaBAGE. LPA long-read matrix (Amersham Biosciences, Piscataway, NJ, USA). 
PGR primers were obtained from MedPrbbe (Oslo, Norway). The primer with GG 
clamp were labeled with 6-carboxyfluorescein (6-FAM) on the tested samples and 

15 tetramethylrhodamine (TMR) on the intemal standards. Primer used for SNP in the 
lymphotoxin alpha (LTA) gene (accession number 153440, locus 6p21.3, NGBI 
reference SNP ID: rs 909253), 5' CGC GGG GGG GGG GGG GGG GCG GTG GGG 
GGG GCG GGG GGG GGT GGT GGG TTT GGT TTT GG 3 ' (SEQ ID NO: l)and 
5' GAG GAG AGG GAG AGA GAG AGA G 3' (SEQ JD NO:2). Primer used for 

20 mutation in Adenomatous polyposis coli (APG) gene (accession nimiber 175 100, 
locus 5q21-q22, Ip34.3-lp32.1, exon 15, codon 1450), 5'- GGG GGG GGG GGG 
GGG GGA GGG GGG GGG GGG GGG GCG GGG GAG GAT TTA CTG GAG 
GTT GGT 3' (SEQ ID NO:3)and 5' AGG TGG TGA AAG AGG TGA AA 3' (SEQ 
ID NO:4). Primer used for SNP in the breast cancer-2 (BRGA2) gene (accession 

25 number 600185, locus 13ql2.3, NCBI reference SNP ID: rs573014), 5' CGC GCG 
GCG CGC GGG GGG GGG GTG CCG GCG GCG GGG CGC GAA GG TAT GTG 
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CAT TGT TTT T 3 ' (SEQ ID N0:5) and 5 ' CCG CAA TAA AGC AAA TAT TAG 
3' (SEQIDNO:6). 



Capillary Electrophoresis (CE) i 

All CE experiments were performed on an MegaBACE 1000 96-capilla3y 
5 DNA analysis system (Amersham Biosciences). The instrument was equipped with 
an automated loading robot Caddy 1000 (Watrex Praha, Prague, Czech Republic), to 
allow for unattended automated operation. To reach high temperatures needed for 
some mutant separation, the temperature sensor was equipped with an additional 
resistor resulting in a positive offset of temperatures by approximately 10 °C. The 
1 0 temperature was recorded using Fluke logging thermometer with Fluke View 

software (Michell Instrument, San Marcos, CA, USA). The temperature profiles 
were constructed using MegaBACE configuration selector (MBCS) software 
(Genomac International, Prague, Czech Republic). The data was processed by SNP 
Profiler, software (Amersham Biosciences). 

15 PGR reaction 

Full blood from blood donors at Ulleval Hospital (Oslo, Norway) was 
anonymously collected and genomic DNA was extracted with QIAamp DNA Blood 
Mini ICit fi*om Qiagen (Valencia, CA, USA). All reactions were performed on a 
PTC-200 thermocycler (MJ Research, Waltham, MA, USA), by mixing 50 ng of 

20 genomic DNA with 25 mM of each dNTP (Abgene, Epsom, UK), 10 Tag buffer, 1 
unit of Tag polymerase (Applied Biosystems, Foster City, CA, USA) and 5 pmol of 
each primer (MedProbe) in a final volume of 25 L. Same cycling conditions were 
applied for amplification of all fragments. The cycling program included 
denaturation for 60 s at 94 °C, annealing for 60 s at 53 °C and elongation 60 s at 72 

25 ^C, for 35 cycles. Where applicable, heteroduplexes were formed by heating the 
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PGR products at 94 for 5 min, then annealing the fragments at 65 °C for 60 min 
followed by slow cooling to 4 °C. 



Optimization of cycling temperature range 

For each target sequence, extended by the GC clamp, the theoretical melting 
5 temperature profile was fibrst calculated using WinMelt. simulation program 
(Medprobe) based on Poland's algorithm. A melt profile will show regions of 
theoretical high and low melting domains of a known sequence. Location of primers 
and GC clamps can be optimized by analyzing their effect on the overall fragment 
melting profile. The temperature oscillation was within ± 1 °C range from the 
10 melting temperature of the low-melting domain (targeit sequence). 

For separation of wild-type and mutant fragments based on differential 
melting, this condition applies. In most of these cases, the time difference between 
the fastest migrating peak from tlie mixture (usually the unincorporated primer) and 
the slowest migrating peak (usually the most denatured fragment) is typically around 
15 10 min relative to the total analysis time of 30 - 40 min (Kristensen, A. et al^ 2001, 
BioTechniqiies, 33:650-654; Zhu, L. et al, 2001, Electrophoresis, 22:3683-3687). 
The multiple-injection scheme could be directly applied to DNA separation, where 
samples would be serially injected and separated at a constant temperature (FIG. 
lA). 

20 A diflFerent situation occurs if a temperature gradient needs to be applied. In 

temperature gradient capillary electrophoresis (TGCE), the running temperature is 
gradually changed during the entire run. The fundamental assumption is that the 
samples have to be subjected to a proximity of optimum melting conditions over a 
sufficient duration of their migration in the capillary. With multiple injections, the 

25 samples eluting at the beginning would be subjected to different temperature ranges 
compared to the ones injected later as shown in FIG. IB. In order to subject samples 
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firom all injections to identical melting conditions, the temperature has to be 
periodically changed to follow the profile of repeatedly injected samples. With this 
arrangement, each sample xmdergoes the same number of temperature gradient 
cycles and thus is subjected to the same melting conditions (FIG. IC). 
5 Others have subjected samples (PGR products) containing an artificial 

high-melting domain ("CG-clamp" or "GC-clamp") to a temperature gradient to 
cover a range of optimum melting conditions (Kristensen, A. et al^ 2001, 
BioTechniques, 33:650-654). A typical result of this single-sweep gradient 
experiment (TGCE) is shown in FIG. 2A. An LTA mutant was subjected to a 

10 descending temperature gradient starting at 52 and ending at 48 °C with a rate of 
0. 1 °C per min. During the gradient, the period of time at which the samples were 
subjected to their optunum melting conditions is given by the gradient slope and has 
direct impact on separation resolution (Minarik, M. et aL, 2001, poster presented at 
ASHG Meeting, San Diego, CA). Described herein are methods for subsidizing this 

15 period with a series of cycles during which the separated samples are several times 
subjected to the melting optimum. FIG. 2B shows a separation of the same sample 
(LTA) subjected to a cycling temperature gradient. It can be seen that the resolution 
in this case is fully comparable to the single-sweep gi*adient experiment from FIG. 
2A. It seems that the overall retardation effect due to partial denaturation in case of 

20 a cychng gradient is similar to a single-sweqp gradient. 

Example 2, Direct identification of all oncogenic mutants in KRAS exon 1 by 
cy^cling temperature capillary electi'ophoresis. 

Over the past few decades, advances in genetics and molecular biology have 
revolutionized the understanding of cancer initiation and progression. Molecular 
25 progression models outlining genetic events have been developed for many solid 
tumors, including colon cancer. Previous reports in the hterature have shown a 
relationship between different KRAS mutations and prognosis and response to 
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medical treatment in colon cancer patients. Furthermore, the presence of a mutated 
KRAS has been correlated with different clinicopathological variables including age 
and gender of patients and tumor location. The mutation analysis method described 
herein is adapted to a 96-capillary electrophoresis instrument that allows 
5 identification of all 12 oncogenic mutations in KliAS exon 1 under denaturing 
conditions. To determine the optimal parameters, a series of DNA constructs 
generated by site-directed mutagenesis was analyzed and the migration times of all 
mutant peaks were measured. A classification tree was then made based on the 
differences in migration time between the mutants and an internal standard. A 
10 randomized series of 500 samples constmcted with mutagenesis as well as 60 blind 
samples firom sporadic colon carcinomas was analyzed to test the method. No 
wild-type samples were scored as mutants and all mutants were correctly identified. 
Post polymerase chain reaction (PGR) analysis time of 96 samples was performed 
within 40 min. 

•\ 

15 DNA samples 

Control mutants and the series of 500 test samples were made by PCR 
mutagenesis. Genomic DNA fi-om the cell line HT 29 with wild-type sequence in 
KRAS exon 1 (ATCC cell repository, ATCC HB 8245) was used as template. One 
sense mutagenesis primer for each missense mutation in codon 12 (wild-type 

20 sequence GGT) and 13 (wild-type sequence GGC) of KRAS exon 1 was designed. 
The 12 mutagenesis primers were 45 base pairs long each and spanned codon 1-15 
of KRAS exon 1. Each primer had a one base difference fi-om the wild-type 
corresponding to the 12 different mutations. The sequence of the antisense primer 
was as follows: 5' TAMRA CGC CCG CCG CGC CCC GCG CCC GTC CCG 

25 CCG CCC CCG CCC GCC TCT ATT GTT GGA TCA TAT TC 3 ' (SEQ ID NO:7; 
products used for internal standard) or 5' 6-FAM CGC CCG CCG CGC CCC GCG 
CCC GTC CCG CCG CCC CCG CCC GCC TCT ATT GTT GGA TCA TAT TC 
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3' (SEQ ID NO:7; control samples and series of 500 test samples). All mutagenesis 
PGR products were sequenced on a MegaBACE 1000 (Amersham Pharmacia 
Biotech, Oslo, Norway) to verify correct sequence, with conditions recommended by 
the manufacturer and 5' CGC CCG CCG CGC CCC GCG CC 3' (SEQ ID NO:8) as 
5 sequence primer. DNA from 60 sporadic colon cancers was extracted with a 
standard DNA extraction kit (Qiagen, Valencia, CA, USA). 

Thermodynamic 

KRAS exon 1 with attached GC-clamp was analyzed with computer program 
MacMelt (MedProbe, Oslo, Norway). The program calculates the thermodynamics 
10 of the transition of double-stranded to single-stranded DNA (Fixman, M. and Freire, * 
J., 1977, Biopolymers, 16:2693-704) 

PGR 

All PGR reactions were performed on a PTC-200 thermocycler (MJ 
Research, Waltham, IMA, USA). The cycling parameters were 35 cycles of 

15 denaturation for 1 min at 94 anneahng at 53 ""C for 1 min and elongation at 72 
for 1 min, followed by heteroduplex fomiation by boiling and re-annealing at 65 
°C for 30 min. PGR amplification was performed by mixing 50 ng genomic DNA or 
0.01 iLiL of PGR product with 25 pM of each dNTP (Perkin Ehner, Oslo, Norway), 
10 Tag buffer, 1 unit cloned Taq (Stratagene, La Jolla, CA, USA), and 5 pmol of 

20 each primer (MedProbe) in a final volume of 25 \xL. PGR primers for amplification 
of KRAS exon 1 in the tumor samples were as follows: 5' ATG AGT GAA TAT 
AAA GTT GTG 3' (SEQ ID NO:9) and 5' 6-FAM GGC CCG CCG CGC CCC 
GCG CCC GTC CCG CCG CCC CCG CCC GCC TCT ATT GTT GGA TCA TAT 
TC3'(SEQIDNO:7). 



wo 2005/047541 PCT/IB2003/005304 



-26- 

Standard 

An internal standard with five peaks was created by mixing TAMRA-labeled 
PGR products Jfrom three homozygous mutant samples (a codon 12 GGT to TGT 
mutant, a codon 12 GGT to CGT mutant and a codon 13 GGC to AGC mutant) with 
5 the unlabeled wild-type PGR product. The product mixture was thereafter boiled 
and re-annealed for heteroduplex formation. Since only the mutant samples were 
labeled prior to mixing and two heteroduplexes co-migrated, five peaks were visible 
when the standard was analyzed with GTGE. The standard wiEis added to all wells in 
the plates prior to GTGE analysis of test samples. 

10 GTGE 

Samples were analyzed by GTGE in a standard 96-capillary DNA sequencing 
instrument, the MegaBAGE. 1000 DNA Analysis System (Amersham Biosciences, 
Uppsala, Sweden). The distance firom the anode to the detector was 40 cm. Regular 
3% linear polyacrylamide (LP A) with urea (7 M) was replaced prior to every run by 

1 5 applying high-pressure nitrogen. PGR products were diluted 1 :25 in water prior to 
electrokinetic injection, accomplished by applying 10 kV for 12 s. The 
electrophoresis condition was a constant field of 150 V/cm. The partial melting of 
DNA firagments was achieved by a combination of urea in the gel and temperature 
surrounding the capillaries. The separating temperature used was a gradient going 

20 fi-om 51 °G down to 49 °G in 1 "^C increments. Each temperature was held for 30 s. 
The gradient was cycled 15 times before the product reached the detector. The 
gradient was followed by constant temperature at 44 °G until all peaks had eluted. 
All samples were analj'zed in 96-well plates. In the presented case, laser-induced 
fluorescence detection was used with excitation at 488 nm (blue laser) and detection 

25 of emission firom the 6-FAM-labeled PGR primer at 520 and excitation at 532 nm 
(green laser) for the TAMRA-labeled PGR primer with detection of emission at 580 
nm. By using these PGR primer labels no spectral overlap was observed between 
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the channels. Migration times of all peaks ia the electropherograms were scored by 
computer program Genetic Profiler and the relative migration times between peaks 
were calculated. 

Classification tree 

5 A plate with 96 control samples (12 mutants times 8) was used to inake a 

classification tree. Each well in the plate contained the TAMRA-labeled standard 
and one 6-FAM-labeled control sample. The plate was analyzed by CTCE under 
various denaturing conditions where temperature and nimiber of cycles were 
adjusted imtil optimal reproducibility was obtained. At denaturing conditions with 

10 high reproducibility, the migration times of wild-type and mutant peaks, as well as 
migration times of the five standard peaks were measured in each electropherogram 
by Genetic Profiler computer software (Amersham Biosciences). Migration times 
fi-om eiglit runs for each mutant were used as references. The migration time data 
was used to generate a classification tree by extrapolating firom the compiled 

15 empirical data (FIG. 9). Identification of the unknown mutants was achieved by 

comparison of migration differences between sample and standard peaks and related 
to the parameters in the classification tree. It is crucial that the classification tree is 
followed fi'om top to bottom for correct identification of the mutants. The letters in 
file classification tree boxes indicate different peaks as noted in FIG. 9. For 

20 example, a GGT to TGT mutation in codon 12 of KRAS has the following 
characteristics: (i) the first heteroduplex (HI) has a shorter elution time than 
standard peak 4 (S4), (ii) the migration time of S4 minus migration time of HI is 
less than 12, and (iii) the migration time for the mutant 

pealc (M) is longer than the migration time for the wild-type peak (W). The flow 
25 chart was used to identify the mutant sequence in 500 blindly scored samples and 60 
sporadic colon carcinomas with knov^m mutation status. 
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Results and discussion 

A theoretical melting profile analysis of the 151 base pair target fi-agment, 
including exon 1 of KRAS and the PCR-attached GC-clamp was first performed. 
All oncogenic mutations in KRAS exon 1 as well as wild-type were analyzed. The 
5 relative differences in thermodynamic properties between the mutations and 
wild-type fragments are illustrated in FIG. 7. Wild-type KRAS exon 1 with the 
sequence GOT and GGC in codons 12 and 13, respectively, is coded as 100% m 
FIG. 7. All but one G to C mutation proved to be more stabile than the wild-type 
sequence, and thus separated out before wild-type peaks in electropherograms. A 
10 GGC to GCC mutation, which was only slightly more stable than wild-type, co- 
migrated with the wild-type. All other mutations were less stable than wild-type and 
eluted after the wild-type peak. Some of the mutations (e.g., AGO and AGT, and 
TGT, TGC and GTT) had very similar melting profiles. 

Based on the theoretical thermodynamics of the different firagments a 
1 5 standard was designed. Three mutations with similar properties to other mutants 
were used in the mtemal standard. The inclusion of these mutations in the standard 
allowed for identification of mutants with almost equivalent themaodynamics. The 
second step required for the CTCE setup was tlie determination of the optimum 
denaturing conditions for the separation of homoduplexes and heteroduplexes by 
20 performing several runs with different numbers of temperature cycles and different 
temperature intervals. The goal was to find denaturing conditions that allowed 
separation corresponding to the theoretical thermodynamics of all mutants. Highly 
reproducible separation of mutant homoduplex peaks was achieved by a gradient 
ranging fi-om 51 down to 49 °C in steps of 1 *=*C. Each temperature was held for 
25 30 s. The gradient was cycled 15 tunes before the product reached the detector. 

Due to local temperature differences in the capillary chamber, few capillaries will be 
exposed to the optimal separating temperature if the temperature was held constant 
during electrophoresis. By cycling the temperature aroimd the theoretical optimal 
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separation temperature, all capillaries will be exposed to the "optimal" temperature, 
during the temperature gradient created in the capillary chamber during 
electrophoresis. The more times the gradient is cycled the longer each capillary will 
be under optimal separation conditions. If the gradient is cycled to too rapidly the 
5 heat capacity of the chamber will not have time to adjust, thus making suboptimal 
denaturing conditions. By cycling the gradient few times, or keeping the 
temperature at the same interval in the cycle for a longer time, suboptimal 
denaturing conditions are observed in some capillaries. Methods for CDCE of 
KRAS exon 1 appUed to a MegaBace 1000 had been reported (Bj0rheim, J. et al, 

10 2002, Anal Biochem,, 304:200-205). Although the methods proved sensitive and 
specific in mutation detection it was not possible to directly identify the mutations 
by this method. Consequently, mutation verification had to be performed by mixing 
of samples with controls prior to re-analysis. 

With the selected denaturing conditions eight 96-well plates with standard 

15 and different mutations and wild-type samples were analyzed by CTCE. Separation 
between the wild-type homoduplex and the heterozygous peaks was achieved in all 
mutated samples. Separation between the wild-type homoduplex and mutant 
homoduplex was found in 1 1 out of 12 mutants. The mutant homozygous peak of 
GGC to GCC transversion in codon 13 co-migrated with wild-type at these 

20 denaturing conditions. These observations corresponded with the theoretical 
thermodynamics of the different fragments. Electropherograms of all mutant 
samples and the standard are shown in FIG. 8. 

No spectral overlap was observed between the channels. Individual 
mutations proved to have distinct migration times when compared to the peaks of 

25 the internal standard and other mutations. Migration times that were able to separate 
different mutations were selected and a classification tree was made (FIG. 9). For 
example, a TGT mutation was characterized first by a shorter migration time of HI 
than S4, and migration time of S4minus migration time of HI was less than 12. 
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Finally, the migration time of the mutant homoduplex (M) was longer than for the 
wild-type peak (WT). No other mutations had these characteristics. 

To test whether blind samples could be scored and identified based on peak 
migration tunes and the classification tree, 500 samples with wild-type and mutant 
5 sequence in codon 12 and 13 of KRAS were scored. The samples were PGR 
amplified by another investigator and analyzed without any knowledge of the 
samples. All 500 samples analyzed were scored correctly. No wild-type samples 
were scored as mutants and no mutant sample was placed in the wrong group. To 
further investigate the adaptability of the method to clinical use, 60 sporadic colon 

10 carcinomas were analyzed. A total of 22 samples were found mutated in KRAS 
exon 1 and all mutants were scored correctly when compared to previous CDCE 
analysis of the samples. Among the 560 samples used to evaluate the method, 100% 
sensitivity and specificity were demonstrated. 

Other reports on mutation analysis of KRAS with automated CDCE have 

15 demonstrated a detection limit as low as 0.1% for heterozygous peaks and down to 
1% for homozygous mutant peaks. Similar detection limits were found in this setup 
(data not shown). Smce tlie decision maldng tree is dependent on mutant 
homoduplex peaks for mutation identification, mutants with mutated allele fraction 
below 1% are not identified in this system. To identify a mutant below tlie 1-% 

20 level, the mutant was mixed with a control mutant, denatured and thereafter 

re-annealed prior to CTCE analysis (Guldberg, P. and Guttler, F., 1993, Nucleic 
Acids Res,, 21 :226 1-2262) as reported for other melting gel techniques. Evaluation 
of the detection limit was performed by the formula as described by Ekstrom et al 
(2000, BioTechniqties, 29:582-589). 

25 Many metliods have been employed to detect mutations. In one approach, 

PGR is part of the detection system. Two widely used methods, allele-specific 
amplification and mutant-enriched PGR, can detect one mutant allele against a 
background of 10 000 wild-type alleles (NoUau, P. and Wagener, C, 1997, Clin, 
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Chem,, 43:1114-1128; Andersen, S. etaL, 1999, Gut, 45:686-692.). Although the 
sensitivity of such methods is high, they are not amenable to automation and high 
throughput. Li a second type of approach, mutations are analyzed after the target 
sequence is amplified by PGR and detected using techniques such as ligation. The 
5 main disadvantage of ligation is that three pairs of oligonucleotides are needed to 
identify a bi-allelic marker (Lehman, T. et al, 1996, Anal Biochem., 239:153-159.). 
Interestingly, PNA-directed PGR in combination with MALDITOF reports high 
sensitivity in KRAS detection, and the possibility for automation. However, several 
steps are needed to identify mutant samples and identification of a possible mutant 

10 in codon 12 must be evaluated prior to codon 13-mutation status analysis (Sim, X. et 
aL, 2002, Nat. BiotechnoL, 20:186-189). ' 

The MegaBAGE 1000 with 96 capillaries in parallel analyzed 96 samples 
under denaturing conditions within 40 min. A very limited number of steps were 
needed m this setup. The samples were PGR-ampUfied in 96-well plates and diluted 

15 1 :25 in water. Thereafter, without any modifications, the plates were ready for 

GTGE analysis. All codon 12 and 13 mutations were directly scored and identified, 
thus eliminating the need for further evaluation. 

GTGE methods of the present invention offer two major advantages over 
previously reported techniques adapted to KRAS mutation analysis. The first 

20 advantage allows for direct identification of all oncogenic KRAS mutations in exon 
1 without the need for sequencing. Secondly, the short post-PGR analysis time 
makes this technique highly desirable for fully automated mutant screening. 

Example 3, Population Screening of Single-Niicleotide Polymorphisms Exemplified 
by Analysis of 8000 Alleles. 
25 Described herein is a method in which the population firequency of 

single-nucleotide polymorplaisms (SNPs) can be efficiently detected and their allele 
firequencies accurately measured. Selected SNPs in TNFp, IL-4, and GTLA-4 were 
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used to demonstrate the method. Blood from 4000 individuals was pooled, DNA 
was extracted, and target sequences were PCR-amplified and anal^'zed by denaturant 
capillary electrophoresis. Alleles were separated into peaks based on melting 
properties of the double DNA helix. Frequencies of the different alleles were 
5 determined by calculating the area under the peaks. Allele frequencies and 
Hardy- Weinberg equilibrium estimated from the pooled data were verified by 
analyzing 7.5% of the samples randomly selected from the blood donor series. The 
method herein is suitable for single-samples and/or pooled-samples analysis of 
SNPs, in which sample treatment is kept to a minimum. The potential throughput of 

10 the method is beyond obtainable numbers of samples. 

SNPs are considered to be useful polymorphic markers for genetic studies of 
pharmacogenetics and polygenic traits (McCarthy, J. and Hilfiker R., 2000, Nat 
Biotechnol, 18:505-508; Roses, A., 2000, Nature, 405:857-865). A worldwide 
effort to collect SNPs has achieved an accumulation of millions in public databases. 

15 However, most of these SNPs have been identified by examination of a limited 
number of individuals, and information on their allele frequencies is lacking or 
tentative. Furthermore, studies have shown that allele frequencies vary widely 
between different ethnic populations. Thus, validation of SNPs and estimation of 
their allele frequencies, especially for each etlinic group, are required before these 

20 markers can be used for genetic studies (Risch, N., 2000, Nature, 405:847-856). 

Pooling of samples is an obvious approach for frequency determination of 
SNP in large populations because it drastically reduces the cost of the analysis and 
labor time compared with genotyping individuals and coimting alleles. The 
quantification method in the pooled analysis must accurately reproduce the genotype 

25 frequencies found when analysis of single samples is performed. 

Described herein is a method for allele frequency estimation of SNPs, based 
on, for example, in part on denaturant capillary electrophoresis (DCE). DCE is 
based on the melting gel theory described by Fischer and Lerman (Fischer, S. and 
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Lerman L., 1983, Proc. Natl Acad. Sci. USA, 80:1579-1583). Double-stranded 
DNA firagments melt characteristically based on the nucleotide sequence and length 
of the fragment when exposed to denaturants such as temperature, urea, or 
formamide. Fragments with only minor differences to the wild-type, such as a base 
5 substitution, result in different melting characteristics. Due to the high sensitivity 
and simpUcity described for melting gel techniques applied on capillary 
electrophoresis (Khrapko, K. et al, 1994, Nucleic Acids Res,, 22:364-369) this 
method was found suitable for population frequency analysis of SNP. 

Whole blood from 4000 donors was prepared in 4 pools, with 1000 samples 

10 in each pool, prior to DNA extraction. The pooled samples and 7.5% of the 

individual DNA samples were subjected to PGR subsequently followed by DCE 
analysis. Allele frequencies were determined in the pooled samples by measuring 
the area under the peaks in the electropherograms. The individual samples were 
genotyped by comparing the peak pattern with a heterozygous standard run 

15 simultaneously with a different fluorochrome. The throughput of this SNP approach 
is limited by the number of samples obtainable. In simple terms, 96 pools with 1000 
samples in each, totaling 96,000 samples, can be analyzed in 35 minutes. For 
single-samples analysis, the theoretical throughput is about 46,000 samples in 24 
hours in a ftilly automated system. This simple 2-step protocol, PGR and 

20 electrophoresis on standard capillary DNA sequencing equipment, could make a 
significant contribution to SNP analysis. 

Pooled samples 

Blood samples from 4000 blood donors were collected from August to 
October 2000 within a 10-week interval to exclude parallels from the same 
25 individual. The sampling was anonymous; no record exists that can link the blood 
sample to the donors. White blood cells were not coimted due to the statistical 
observation described in the central limit theorem. 
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Z- ' is fipjiTOixiiiiateJy NfOJX 

is approximately N(0,1). 

The theorem states that in a random samphng of a population with mean ii 
and standard deviation of a, the distribution of average X when n is large is 
approximately normal, with mean \i and standard deviation c/square root of n 
5 (Bhattacharyya, G. and Johnson. R., "The normal distribution and random samples." 
In Bradley RA, Kendall DG, Hunter JS,Watson GS (eds): Statistical Concepts and 
Methods. New York: John Wiley, 1977:187-232). Consequently, the average 
number of white blood cells for each allele will approach the average. 

Tw:o pooling steps made 4 pools with blood from 1000 donors in each. First, 

10 50 jiL blood from 100 donors was pipetted into a tube and mixed vigorously. 

Subsequently, 1 mL blood from 10 pools made up of 100 donors was mixed, giving 
4 pools witli 1000 samples in each. DNA was extracted from the pools by use of the 
Qiagen Max Blood extraction kit (Qiagen Inc., Valencia, CA). All single samples 
subjected to DNA extraction from the blood donor samples were extracted with a 

15 GenoJVF'^-48 (Geno vision, Oslo, Norway), and aliquots were transferred into 
96-well microliter plates. 

PGR reactions 

All reactions were perfomied on a PTC-200 themial cycler (MJ 
Research, Waltham, MA) by mixing 50 ng of pooled- or individual- 
20 sample DNA with 25 mM of each dNTP (Perkin Ekiier, Oslo, Norway), lOX Tag 
buffer, 1 unit of Tag polymerase (Perkin Elmer), and 5 pmol of each primer 
(MedProbe, Oslo, Norway) in a final voliune of 25 \iL. The same cycling 
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parameters were used for the 3 fragments- denaturation 60 seconds at 94 
annealing 60 seconds at 53 elongation 60 seconds at 72 *^C- for 35 cycles. 

The PGR products were denatured 5 minutes at 94 "^C and incubated 30 
minutes at 65 °C for heteroduplex formation. 

5 Primers 

TNFP, Genbank accession number M55913, variation in base pair 329 A/G: 
5' CGC CCG CCG CGC CCC GCG CCC GTC CCG CCG CCC CCG CCC GCT 
GGT GGG TTT GGT TTT GG 3' (SEQ ID NO:l); 5' GAG GAG AGG GAG AGA 
GAG AGA G 3' (SEQ ID NO:2). 
10 IL-4, Genbank accession number AF395OOS5 variation iii base pair 12559 

TIC: 5* CGC CCG CCG CGC CCC GCG CCC GTC CCG CCG CCC CCG CCC 
GCT ATC TTT GTC AGG ATT GGA T 3' (SEQ ID NO:10); 5' ATG CTA GCA 
GGAAGAACAGA3' 

GTLA-4,' Genbank accession number AF414120, variation in base pair 204 
' 15 A/G: 5' CGC GCG CCG CGC CGC GGG CGC GTC GGG GCG CCG GCG CGC 
GGT TGC TGA AGA GGT GAA GAG 3' (SEQ ID NO:l 1); 5' GAG GGA TGA 
AGA GAA GAA AA 3' (SEQ ID NO: 12). 

The primers with GC-clamp were labeled with 6-carboxyfluorescein or 
Tamra fluorochromes. 

20 Denaturant capillary electrophoresis 

The pooled samples were analyzed by the ABI 310 Genetic Analyzer 

(Applied Biosystems, Foster City, CA) under constant denaturing conditions. 

Temperatm-e of the capillary was set to 52 °G, 47 °G, and 56 °G (TNFp, rL-4, 

CTLA-4, respectively). 
25 Single samples were analyzed by temperature gradient capillary 

electrophoresis (TGCB) in a 96-capillary system. The instrumentation was a 
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staiidard multi-capillary DNA sequencing instrument MegaBACE'^ 1000 DNA 
Analysis System (Amersham Pharmacia Biotech, Oslo, Norway). The instrument is 
a high-throughput, fluorescence-based DNA system using capillary electrophoresis 
witti 96 capillaries in parallel. The distance jfrom the anode to the detector is 40 cm. 
5 The MegaBACE™ platform uses cofocal laser scanning focused on each capillary. 
Regular 3% linear polyacrylamide (LP A) with urea was replaced prior to every run 
by applying high-pressure nitrogen. PCR products were diluted 1 :25 in water prior 
to electrokdnetic injection, accomplished by applying 10 kV for 12 s. The 
electrophoretic condition was a constant field of 150 V/cm. The partial melting of 

1 0 DNA fragments was achieved by a combination of urea in the matrix and 
temperature surrounding the capillaries, 

TNpp was separated with a gradient going from 54 °C down to 52 in 
steps of 0.5 *^C, and each temperature was held for 210 seconds. The gradient was 
followed by constant temperature at 50 ®C imtil the end of the run. 

15 The gradient used to separate the C and T alleles in IL-4 was 51 *^C, 50 "^C, 

and 49 °C, and each temperature was held for 2 minutes. The gradient was cycled 4 
times, followed by constant temperature at 48 °C imtil the end of the iTin. 

CTLA-4 alleles were separated by a cycling gradient of 56 ^C, 55 ""C, and 54 
(1 min each temperature) 7 times, followed by constant temperature at 50 ''C. 

20 The areas under the peaks were measured with the Aclaiowledge computer 

program (Biopac Systems, Inc., Santa Barbara, CA). To detect potential PCR bias 
between the 2 alleles, peak areas in individual heterozygous samples were measured. 

Selected samples were sequenced to verify the genotype, using the first 20 
bases of the GC-clamp primer. The sequencing reaction and analysis were set up as 

25 described by the manufacturer (Amersham Pharmacia Biotech, Oslo, Norway). 
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Control of PGR bias 

Analysis of PCR-amplified DNA introduces the possibility of preferential 
amplification of different alleles. PGR bias was analyzed by measuring the peak 
areas in heterozygous samples. No preferential amplification of alleles was 
5 observed when areas under the peaks were measured in heterozygous individual 
samples. FIG. 1 1 represents 5 individual heterozygous samples in TNFp analyzed 
by TGGE in different runs; the average ratio between the peak area was 0.996 (one 
SD = 0.028). The average ratio between the 2 alleles in IL-4 and GTLA-4 was 0.994 
(one SD = 0.06) and 0.993 (one SD = 0.06), respectively. 

10 Pooled Samples 

In the four pools, each allele jfrequency was detemiined by measuring the 
pealc areas representing each allele. FIG. 2 shows the electropherograms of the four 
pools analyzed for the selected TNPP polymorpliism. The results firom the pooled 
blood samples are sxmmiarized in FIG. 12. For the TNFp polymoiphism, the 

15 frequency was 63. 1% and 36.9% for A alleles and G alleles, respectively. The mean 
allele frequency for IL-4 was 84.6% for C alleles and 15.4% for T alleles. For 
CTLA-4, the frequencies found were 55.1% for A alleles and 44.9% for G alleles. 

Individual samples 

Genotypes of 300 individual samples were determined by coelution with an 

20 internal heterozygous carboxytetramethylrhodamine (Tamra)-labeled standard. In 
FIG. 13, the genotypes in the TNFp SNP are depicted with GG at the top, AA in the 
middle, and GA at bottom. The peaks presented by the bold line are a 
Tamra-labeled standard run simultaneously in each capiUa^5^ There was no need for 
spectral overlap corrections due to the use of 488-nm and 532-mn laser excitation of 

25 6-caroxyfluorescein and Tamra, respectively. In FIG. 12, the allele frequencies for 
the 3 SNP measured in the pooled samples are compared with the single-samples 
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analysis. For TNFp polymorphism, frequency was 63.7% and 36.3% for A alleles 
and G alleles, respectively. The mean allele frequency for IL-4 was 84.8% for C 
alleles and 15.2% for T alleles. For CTLA-4, the frequencies found were 55.7% for 
A alleles and 44.3% for G alleles. The Hardy-Weinberg equiUbrixun was calculated 
5 from the allele fi*equencies obtained from the pooled data. The estimated genotype 
distributions were foxmd -equal to the observed genotype frequencies from the 
single-samples data. 

Sample throughput 

Multiple injection of samples during electrophoresis increased tlie 
10 throughput of the method. New samples were injected in each capillary with a 8.5 

minute interval. FIG. 14 displays the electropherogram after injections of 5 samples. 
To test whether the method was amenable on a 384-capillary instrument, the 

samples were re-analyzed on a MegaBACE™ 4000. The denaturing conditions 

were controlled with the same gradient used on the 96-capillary instrument. No 
15 apparent differences were observed between the instruments (data not shown). This 

opens a theoretical throughput of 7680 samples per hour (5 injections x 4 channels x 

384 wells). 

Discussion 

SNP discovery has resulted in large databases containing information on 
20 DNA variation in miUions of sequence positions of the human genome. Most SNPs 
reported are detected by analyzing a limited number of individuals, and the 
distribution frequency of an SNP is often based on small populations. Furthermore, 
it is likely that a broad part of the SNPs reported has no or limited impact within 
medicine and biology. SNPs that turn out to be correlated to aberrant observations 
25 are described in several reports and most likely will be added to the list of important 
SNPs in years to come. Medical intervention could be dependent on SNP 
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knowledge so that severe side effects and inefficient treatment can be avoided. 
Fvjither knowledge of SNP allele Jfrequency distribution in different populations will 
be needed for optimal treatment of the individual patients and for certain subgroups 
of patients. 

5 Described herein is a method for population fi-equency analysis of SNP. 

Pooled samples were subjected to denaturant capillary electrophoresis in order to 
measure allele fi-equencies of defined SNPs. Single-samples analysis was performed 
to confirm the allele fi*equencies found in the pooled samples. By pooling large 
groups of samples, several technical problems have to be addressed. Preferential 

10 amplification of one of the alleles during PGR would give a false relative ratio 
between the two alleles. No difference in amplification of the different alleles 
within the target sequences of heterozygous single samples could be observed (FIG. 
10). Ratio of peak areas was close to 1, and thus no preferential ampUfication took 
place. Consequently, peak areas measured in the pooled-samples electropherograms 

15 represent allele copies entered in the PGR reaction (FIG. 11). 

Second, when quantifying alleles by laser-induced fluorescence and 
electrophoresis, the peaks must be within the dj^amic range of the detector, and the 
fluorescence signal should reflect the total copy number. Fragments that differ by a 
base substitution or mismatch (heteroduplexes) will pass the detector with the same 

20 velocity when analyzed below optimal separation conditions (Ekstrem, P. et al^ 

2000, BioTechniques, 29:582-589). Thus, the fluorescence signal will represent the 
number of DNA copies. 

Reports on the analysis of pooled DNA have been restricted to analysis of 
pools made up of equal molar of DNA. By pooling 1000 samples of whole blood or 

25 DNA, no prior Icnowledge of white blood cell co\mt or DNA concentration of the 

individual samples is needed. Given that variance of white blood cell coxmt or DNA 
concentrations is finite, the average number of each allele will approach the mean as 
given by the central limit theorem. 
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The allele frequencies estimated from the pooled data are given in FIG. 12. 
The standard deviation (given in parentheses) reflects the reproducibility of the 
method and the difference between the estimate from each pool. Comparing allele 
or genotype data estimated in the pools with frequencies obtained through 
5 single-samples analysis validates this pooling strategy. By using a standard DNA 
sequencing instrument, information on SNP population frequencies can be obtained 
with niinimum effort. Either 96 different ethnic groups or disease cohorts of interest 
made up of 1000 samples each can be analyzed within 35 minutes. Furthermore, 
this approach was performed on standard sequencing equipment accessible for most 

10 of the research community. The samples processing is simple, involving only the 
almost mandatory PGR step followed by denaturant capillary electrophoresis. There 
is no need for post-PGR handling such as purification, hybridization, or restriction 
enzyme cutting, so hands-on time is taken to a minimum. As a consequence of the 
simple protocol, there is no need for expensive consimiables as reported for other 

15 techniques. As long as the target sequence consists of a high and low melting 

domain (naturally or constructed) and can.be amplified by PGR, virtually any DNA 
sequence can be analj/zed for SNPs with this method. The sensitivity of the method 
has been described previously and opens the possibility of detecting SNPs of all 
frequencies (an SNP is defined to be present in more than 1% of the population) in 

20 pooled samples. 

Single-samples analysis was performed by DGE in a 96-capillary sequencing 
instrument. Initially, constant denaturing conditions were attempted, but large 
capillary variability with respect to allele separation was observed (data not shown). 
Thus, appropriate cycling gradients were used to compensate for temperature 

25 differences within the analysis camber and between the capillaries. Alleles were 

detemiined b}^ coelution with a heterozygous standard (FIG. 13). The clear baseline 
separations of the two alleles in FIG. 13 enable automated peak recognition and 
allele scoring. Single-samples throughput can further be increased by repeated 
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injection of samples iii the same capillary (FIG. 14). Theoretical foundations for this 
approach are described elsewhere (manuscript in preparation). In the laboratory, 
only two main steps have to be addressed: a conventional PGR "amplification of the 
target sequence and DCE analysis. Both steps are ready to be automated by 
5 liquid-handling robotics and automation of capillary DNA sequencing instruments. 
The Hmited number of steps needed for analysis reduces the analysis time, and the 
risk of introducing false results is lowered. 

No new method introduced in SNP analysis can omit the question about the 
potential cost in the case of high-throughput analysis. The pooling strategy reduces 

10 the cost 1000-fold as compared to single-samples analysis. Furthermore, the method 
presented herein can reduce the cost of single-samples PGR reaction 100-fold. The 
DNA sequencing instrument needs to handle 5- to lO-jiL volumes; when taking into 
account the 25-fold dilution factor in HjO of the PGR product, it would be sufficient 
for analysis purposes only to use the 200-nL PGR product. Gonsequently, PGR cost 

15 can be reduced significantly with the use of nano PGR equipment. Furthermore, the 
electrophoresis cost can be reduced 20-fold by using the instrument to its fiill 
capacity by multiple injections and using all four channels. Clearly, this simple and 
inexpensive protocol may be advantageous in analyzing large cohorts. 

Example 4. DNA variants in the ATM gene are not associated with sporadic rectal 

20 cancer in a Noi-wegian population-based study. 

A large nmnber of DNA single-nucleotide polymorphisms (SNPs) have been 
discovered following the Human Genome Project. Several projects have been 
launched to find associations between SNPs and various disease cohorts. Described 
herein is an examination of the possible association between the reported SNPs and 

25 sporadic rectal cancer. It has been proposed that SNPs in the ataxi-telangiectasia 
mutated (ATM) gene modulate the penetrance of some cancers. The investigated 
target sequence harbors three polymorphisms (IVS38-8 T/G in intron 38, 5557 G/A 
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and 5558 A/T in exon 39), resulting in eight possible microhaplotypes at the DNA 
level. Furthermore, the two exonic SNPs are sited next to each other, allowing foxxr 
possible amino acids in the same codon. 

Using CTCE as described herein, SNPs and microhaplotypes based on 
5 theoretical thermodynamics and migration of variant fragments are employed. 

Fluorophore-labeled PGR products were analyzed without any post-PCR steps on a 
standard 96 capillary-sequencing instrument under denaturing conditions. More 
than 7000 alleles were microhaplotyped based on peak migration pattems of 
individual samples and sequencing results. The ATM polymorphisms and 

10 microhaplotypes examined did not significantly differ between sporadic rectal 
cancer and normal population. 

The SNP analysis method used in this study employs the following: (a) PGR 
amplification of the target sequence followed by (b) allelic discrimination by DNA 
melting equilibrium. In fliis study, three polymorphisms reported in intron 38 and 

1 5 exon 39 of the ATM gene within a 1 00-bp fi-agment for SNPs and their 

microhaplotypes are analj^ed. A total of 3,526 samples from a normal population 
and 151 samples from sporadic rectal cancer patients were examined to test the 
method and estabUsh microhaplotype frequencies in the respective populations. 

Study population 

20 Tissue samples from 151 patients (90 men, 61 women) with sporadic rectal 

cancer were collected at the Norwegian Radium Hospital between 2000 and 2002. 
Medical records verified a diagnosis of rectal cancer (n =151). Patients' ages 
ranged from 26 to 87 years with a mean age of 65. All samples were collected after 
infomaed consent. None of the patients had a family history of colorectal cancer. 

25 The control population consisted of 3,525 anonymous blood donors serving the 
blood bank at Ulleval hospital, Oslo, Norway. The age distribution and gender are 



wo 2005/047541 



PCT/1B2003/005304 



-43- 

therefore xinkaown. Whole blood samples were collected witiun a narrow time 
frame (10 weeks) to exclude sampling of duplicates. 

DNA extraction 

Tissue was preserved in RNA-later (Ambicon, Austin, Tex., USA) xmtil 
5 DNA extraction was performed. DNA was extracted from tumor tissue with 

QIAamp DNA Kit (Qiagen, Valencia, Calif., USA). DNA was extracted from blood 
by used of GenoM-48 (Geno Vision, Oslo, Norway) magnetic bead DNA extractor. 
Aliquots (100 \iL) of extracted DNA were dispensed into 96-wells plates (Abgene, 
Epsom, United Kingdom) prior to PGR. 

1 0 Polymerase chain reaction 

Primers for the target sequence in ATM gene (GenBank accession number 
U82828.1, bp position 92313.92461) were designed by the primer design program 
Primer 3 (Whitehead Institute, Cambridge, MA), and the thermodynamics of the 
fragments was analyzed by the WhiMelt computer program (Lerman, L. and 

15 Silverstein, K., 1987, Methods Enzymol, 155:482-501). A 60-mer 

carboxyfluorescein (6-FAM)-labeled primer with GC-clamp (5' 6-FAM-GCG GGC 
GGC GCG GGG GGC GGG GAG GGC GGC GGG GGC GGG C-TC AAA CTA 
TTG GGT GGA TTT G 3'; SEQ ID NO: 13) and a 20-mer primer (5' TCC CTG 
AAC ATG TGT AGA AAG C 3'; SEQ ID NO: 14) were used for PGR ampHfication. 

20 A TAMRA-labeled primer (5* TAMRA - GCG GGC GGC GCG GGG CGC GG 3*; 
SEQ ID NO: 15) was used in place of the 5' 6-FAM-labeled primer, to create an 
intemal standard by re-amphfying a heterozygous samples of the two most frequent 
microhaplotypes (IVS38-8T-5557G-5558A, IVS38-8T-5557A-5558A). PGR 
reaction mixtures consisted of 50 ng genomic DNA, 0.06 U/|aL Taq polymerase, IX 

25 Buffer (ABgene) containing 2.5 mM MgClj, and 0.4 mM dNTP mix (ABgene). 
Each of the primers was added to a final concentration of 0.24 mM. The total 
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volume was adjusted to 25 |aL with sterile ion-exchanged MilliQ water (MilHpore, 
Oslo, Norway). Amplification was performed in an air-thermocycler Peltier 
Thermal Cycler PTC 200 (MJ Research, Waltham, Mass., USA), using the 
following cycling conditions: 5 minutes at 94 °C, followed by 35 cycles of 30 
5 seconds at 94 °C, 30 seconds at 55 ''C and 60 seconds at 72 ''C. PCR was 
terminated after alO minutes elongation step at 72 °C. 

Gene nomenclature 

The ATM gene is mapped to chromosome 1 lq22.3, and the target sequence 
used was obtained from GenBank, accession no. U82828.1 . The target sequence 
10 was located within intron 38 and exon 39. The polymorphisms examined 

correspond to NCBI SNP reference numbers rs3092829 (IVS38-8 T/C), rsl801516 
(5557 G/A), andrsl801673 (5558 A/T). 

Cycling temperature capillary electrophoresis 

A standard 96-capillary DNA sequencing instrument MegaBACE 1000 DNA 

15 Analysis System (Amersham Biosciences, Oslo, Norway) was used for cycling 
temperature capillary electrophoresis (CTCE). Coated fused sihca capillaries were 
obtained in sets of 16 capillaries (Amersham). The capillaries have an outer 
diameter of 200 pm and an internal diameter of 75 |Lun. . The distance from the 
cathode to the detector is 40 cm. The MegaBACE platform uses a scanning cofocal 

20 laser for detection of fluorescence in each capillary. Standard sequencing 3% linear 
polyacrylamide with urea was replaced prior to every run, by applying high-pressure 
nitrogen gas. PCR products were diluted 1 :20 in water prior to electrokinetic 
injection, accomplished by applying 200 V/cm for 12 s. Electrophoresis was 
performed at a constant 150 V/cm. The partial melting of DNA fragments was 

25 achieved by a combination of urea in the matrix and decreasing cycling temperature 
going from 48 *=*C to 46 °C in 1 ""C steps. Each temperature was held for 30 seconds, 
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10 



and the gradient was repeated 15 times and followed by a hold at 44 ""C until the end 
of the run. With the present setup 96 samples were analyzed within 40 minutes by 
use of CTCE. 

DNA sequencing 

All samples with aberrant peak pattems compared to the intemal standard 
were sequenced. The imlabeled sequencing primer for tlie forward sequence was a 
primer based on the first 20 bases in the GC-clamp (5* GCG GGC GGC GCG GGG 
CGC GG 3'; SEQ ID NO: 1 5) and the same 20-mer primer used for PGR described 
above was used to generate the confirmatory reverse sequence. Samples were 
sequenced using the MegaB ACE standard protocol and sequencing analysis was 
perfomied on the same instrument under sequencing conditions as described by the 
manufacturer (Amersham). 

Thermodynamics 

All microhaplotype combinations of the three polymorphisms with 
nucleotide number IVS3S-8 T/C, 5557 G/A, and 5558 A/T were simulated using 
DNA melting software WinMelt (MedProbe) based on the Poland algorithm and 
variants thereof. Normalization of microhaplotype melting behavior relative to the 
wild type (e,g., tlie most firequent allele combination) in the nucleotide positions 
IVS38-8T-5557G-5558A was calculated using the following formula: 
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where X and Y are estimated melting temperatures for wild type and microhaplotype 
at given base pair (n), respectively. 

Statistical methods 

The two-tailed test was used to calculate statistically significant (P < 0.05) 
5 diffeences in allele, genotype, and microhaplotype firequencies between the normal 
population and the rectal cancer group. 

Results 

Three ATM polymorphisms, IVS38-8 T/C, 5557 G/A, and 5558 A/T were 
analyzed in 3,526 blood donors and 151 sporadic rectal cancer by CTCE. Alleles 

10 were resolved by cooperative melting of double strand DNA. Partial melting of 
DNA fragments was detected as change in mobility during electrophoresis. FIG. 15 
displays five genotypes of the three SNPs examined. The first three samples (A, B, 
C in FIG. 15) labeled \\dth 6-FAM coeluted with the respective alleles of the 
TAMRA-labeled intemal standard. FIGS. 15D and 15E display two 

15 electrophero grams with alleles that did not coelute with the intemal standard. 
Theoretical melting behavior of each microhaplotype was normalized to the 
wild-type IVS38-S T 5557 G 5558 A (FIG. 16). Microhaplotype IVS38-S T 5557 A 
5558 A has positive value and eluted after the wild-type (IVS38-8 T 5557 G and 
5558 A). The remaimng microhaplotypes have negative values (FIG. 16), which 

20 indicate higher temperature properties as compared to the wild-type. Consequently 
these alleles elute before the wild-type (FIGS. 15D and 15E). Knowledge of melting 
behavior for different microhaplotypes relative to the wild-type in combination with 
DNA sequencing allowed direct determination of microhaplotypes. The 
electropherogram in FIG. 15D shows that one peak coeluted with the IVS38-8 T 

25 5557 A 5558 A peak in the standard and one peak eluting before the standard. 

Sequencing this sample revealed heterozygosity in position 5557 G/A and 5558 A/T 
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(data not shown). This information in combination with nomiahzed melting values 
(FIG, 16) made it possible to identify the microhaplotypes as IVS38-8 T 5557 G 
5558 T and IVS38-8 T 5557 A 5558 A. The electropherogram in FIG. 15E shows 
that both peaks eluted before the standard. This indicates that the sample has a 
5 higher melting temperature than the wild-type. Sequencing verified a triple 
heterozygous sample for the fliree SNPs (data not shown). Consequently the 
microhaplotypes of this sample must be IVS38-8 C 5557 A 5558 A and IVS38-8 T 
5557 G 5558 T. 

The observed genotype frequencies of tlie three examined ATM SNPs 

10 (IVS3S-8 T/C, 5557 G/A and 5558 A/T) and the frequencies of the alleles are shown 
in FIG. 17. The genotypes of the tliree polymorphisms were in Hardy- Weinberg 
equihbrium for both groups (data not shown). 

For SNP IVS38-8 there were no patients foxmd in the sporadic rectal cancer 
group with the genotype combination CC. For SNP 5558 the combinations A/T and 

15 TT were not foxmd in the cancer group, and TT was not foimd in the normal 

population. No statistically significant differences between the two groups were 
found for the SNPs, IVS38-8 T/C (P=0.52), 5557 G/A (P=0.28), and 5558 A/T 
(P=0.51). The allele frequencies in the three SNPs were also statistically 
insignificant (IVS38-8 T/C, P=0.49; 5557 G/A, P=0.78; 5558 A/T, P=0.51). 

20 The frequencies of the microhaplotypes are shown in FIG. 16. The observed 

frequencies of the microhaplotypes were in good correspondence with estimated 
microhaplotypes based on allele frequencies (data not shown). In the 7,052 alleles 
analyzed in the normal population there were three microhaplotypes that were not 
found and two of these microhaplotypes consisted of ATT, in codon 1853, wliich 

25 codes for isoleucine. In the sporadic rectal cancer group with 302 alleles only three 
microhaplotypes were foimd. No statistically significant differences in the frequency 
of microhaplotype were fomid within the two groups (P=0.40). 
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Discussion 

The results described herein, obtained by using a novel method for 
identifying and quantifying polymorphisms and allele frequencies in a population, 
demonstrate no correlation between incidence of rectal cancer and the three ATM 
5 polymorphisms IVS38-8 T/C, 5557 G/A, and 5558 A/T and microhaplotypes. 
Although no association was found between selected SNPs and rectal cancer, the 
method used to analyze SNPs, genotypes, or microhaplotypes exhibits a high 
throughput, sensitivity, and robustness for DNA variation analj^sis. 

Example 5. Evaluatiofi of denaturing conditions in analysis of DNA variants 

10 applied to multi-capillaiy electrophoresis instruments. 

Denaturant slab gel techniques (DGGE and variants thereolQ ^md denaturing 
CE (DCE) have been used for the analysis of mutations and single nucleotide 
polymorphisms by several research groups. DCE applied to commercially available 
capillary electrophoresis instruments has been demonstrated to be sensitive, specific, 

15 and robust for mutation and SNP analysis. However, instruments currently available 
on the market are not designed for DCE in their present forms. In both single and 
multi-capillary instruments inaccurate temperature control units lead to 
irreproducible denaturing conditions, both between runs and from capillary to 
capillary within an electrophoresis run. Cycling a temperature gradient several times 

20 aroimd a theoretically defined optimum for a given DNA fragment results in more 
reproducible and controllable denaturing conditions compared to a constant 
denaturant or single sweep temperature gradient In this example, it is demonstrate 
that wdth the temperature differences observed, the common technique of applying a 
temperature gradient is not sufficient for optimal separation conditions in all 

25 capillaries. Furthermore, it is demonstrated that cycling of the temperature around a 
fragment's theoretically calculated melting temperature during electrophoresis leads 
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to reproducible results in all capillaries. These results clearly show the potential for 
standard commercial capillary instruments in automated cycling temperature CE. 

Constant Denaturant Capillary Electrophoresis (CDCE), introduced in 1994 
(Khrapko, K. et aL, 1994, Nucleic Acids Res., 22:364-369) has been demonstrated to 
5 have high sensitivity as a method in mutation detection. CDCE was performed on 
home-built apparatus and was based on the same separation principles as its slab gel 
analogs, Denaturant Gradient Gel Electrophoresis and Constant Denaturant Gel 
Electrophoresis (Hovig, E. et aL, 1991, Mutat Res., 262:63-71). The separation 
principle utilized by these methods is based on differential melting of DNA (Fischer, 

10 S. and Lerman, L., 1980, Proc. Natl Acad. Set USA, 77:4420-4424; Fischer, S. and 
Lerman, L., 1979, Cell, 16:191-200; Fixman, M. andPreire, J,, 1977, Biopolynters, 
16:2693-2704). hi short, mutated sequences can be separated from wild-type 
sequences on the basis of nucleotide sequence and length. A double-stranded 
fragment undergoes melting into discrete single strand domains when exposed to 

15 denaturing temperature or a chemical denaturant (e. g., urea or formamide). Under 
properly selected partially denaturing conditions, different DNA variants will exhibit 
alternative secondary structures that will allow their electrophoretic separation in a 
sieving electrophoresis matrix. The separation performance can be enhanced further 
if the two strands of the DNA fragment are forced to stay in proximity to one 

20 another during the partial denaturing process. This can be easily achieved by 
extending the fragments at one end by a high-melting domain, for example an 
artificial GC-rich sequence (GC-clamp) (Myers, R. et al. Nucleic Acids Res., 
13:3111-3129). The melting temperature of the clamp is designed to be above the 
optimiom partial-melting temperature, therefore holding the two strands together at 

25 one end. 

Over the past several years, CDCE has been used for detection of various 
low-level frequency mutations. Recently, we have demonstrated the use of 
semiautomated mutation analysis of KRAS using the CDCE technique on 
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commercial single and multi-capillary instruments. In those studies, the separation 
temperature was accurately maintained to achieve high resolution. In the 
multi-capillary CDCE format, a home-made solid-state heater accommodating all 96 
capillaries was necessary to maintain constant separation temperature. With an 
5 unmodified apparatus, the temperature in the capillary chamber could not be 

maintained stable enough to permit uniform denaturing conditions in all capillaries. 

The main purpose of the present study was to demonstrate that neither 
constant temperature, nor siagle temperature gradients, allow optimal temperature 
distribution 

10 across all capillaries ia a 96-capillary format. To eliminate this effect, we 

demonstrated that cycling of the separation temperature around the theoretically 
calculated thermodynamic properties of the relevant DNA fragment during 
electrophoresis allowed reproducible results in all capillaries and over repeated runs. 

DNAandPCR 

15 PGR samples with target sequences in KRAS exon 1 and Lymphotoxin alpha 

promoter were used to demonstrate separation in different capillaries under 
denaturing CE. Primers and PGR conditions for the fragments used have 
been described previously (see Examples 1 and 2). Both the KRAS and 
Lymphotoxin alpha PGR products contain one low melting domain covering the 

20 target sequence and one artificial high melting domain. 

An internal standard with five peaks was created by mixing labeled PGR 
products from three homozygous mutant samples (a codon 12 GOT to TOT mutant, 
a codon 12 GOT to GOT mutant and a codon 13 GGG to AGG mutant) witli 
unlabeled wild-type PGR product. Two of the heteroduplex variants co-migrated in 

25 this standard. 



wo 2005/047541 



PCT/IB2003/005304 



-51- 

Theoretical thermodynamics 

Theoretical melting projSles of the different target fragments were obtained 
from the Poland algorithm and by using the WinMelt computer program, applying 
the Melt87 program (MedProbe, Oslo, Norway). The optimal separation 
5 temperature proposed by tlie WinMelt^oland Web page was adjusted on the basis 
of the concentration of urea in the matrix. For each molar increment of urea 
concentration the temperature was lowered by approximately 3 K. 

CE 

PGR samples were analyzed in a standard 96-capillary DNA analyzer 
1 0 (MegaB ACE™ i OOO DNA Analysis System, Amersham Biosciences, Sunnyvale, 
CA, USA). The distance from the anode to the detector was 40 cm. Linear 
polyacrylamide (MegaBACE LP A) containing 7 M urea was replaced in capillaries 
prior to each run. Standard coated capillaries were used. PCR products were loaded 
into the capillaries by electroldnetic injection at 200 V/cm for 12 s. The 
15 electrophoresis was carried out at a constant field of 150 V/cm. All samples were 
analyzed in 96-well plates. Laser-induced fluorescence was used with excitation at 
488 imi (blue laser) and detection of emission at 520 imi (FAM channel). 

Denaturing conditions and temperature control 

For each run, PCR products from the same target sequence were dispensed to 

20 all wells in a 96-well plate and analyzed xmder various denaturing conditions. The 
denaturing temperature in the capillary chamber was either held constant at 
theoretically calculated optimal value or cycled for varied time intervals around the 
calculated optimal temperature. The denaturing conditions, e.g., the temperature 
profiles, were controlled through the Instrument Control Manager (ICM) software 

25 package (Amersham Biosciences). Temperature sensors were placed in close 

proximity to tlie capillaries. All sensors were placed at the same vertical level, and 
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at different horizontal locations. One temperature sensor was placed at the front part 
of the capillaries (A), one in the middle part (B), and one at the back of the 
capillaries (C). The temperature recording from the three thermocouple probes (K 
type) (Sievert Max, Oslo, Norway) was performed using an analog computer board 
5 (lOtech, Inc., Cleveland, Ohio, USA). The temperature probes were calibrated 
before the temperature recording took place, and the rate of measuremrat was one 
point per second(l Hz). 

Results 

The first step required in a capillary electrophoresis protocol using 

10 differential melting on a defined target fragment, is the determination of the optimal 
melting conditions for the separation of homozygous alleles as given by one of 
several available computer programs utilizing variants of the Poland .algorithm. 
These programs include the Poland algorithm, and the computer software program 
WinMelt/MacMelt designed to calculate thermodynamics of double stranded DNA. 

15 FIG. 1 8 demonstrates the calculated relative stabiUty of two alleles of the 
lymphotoxin alpha fragment at different melting probabilities, and the 
experimentally observed difference in migration {e. g,, scan number variant 1 minus 
scan number variant 2). For the experimental data, the highest temperature not 
causing branching of the double stranded DNA was defined as 0% difference in 

20 stability between the alleles. The lowest temperature with complete branching of the 
low melting domain was defined as 100%. At these temperatures no allele 
separation was expected. Note that the thermodynamics calculated using the Poland 
web page (H) predicted the observed optimal thermodynamic separation (9), as 
compared to the actual behavior of DCE analysis of the fragments. WiiiMelt (F), 

25 however, did not predict the presence of a specific optimal themiodynamic 
separation region between the fragments. Thus, quite large differences in 
thermodynamics were predicted at all melting probabiUties with this program. 
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The thermodynamics (differences for each base pair at 50% melting 
probability) of all different oncogenic mutants in KRAS exon 1 were calculated with 
WinMelt at 50% melting probability. Thereafter, the calculated values were 
correlated with observed migration differences for each mutant. A linear correlation 
5 (r^ = 0.97) was found between migration times of the different variants and the 

theoretical thermodynamics calculated (FIG. 19). The same exercise was performed 
using theoretical melting probabilities obtained from the Poland web page, A linear 
relationship found was j-^ = 0.93. 

To be able to transfer the calculated melting properties of a fragment directly 

10 to CE analysis, the temperature distribution in the capillary chamber was examined. 
Extemal temperature sensors were placed inside the capillary chamber of a 
MegaBACE™ 1000 in close proximity to the capillaries. Several temperature 
profiles were constructed to give the same average denaturing temperature (48 °C) in 
the capillary chamber. The results showing local temperature differences within the 

1 5 capillary chamber in these experiments are presented in FIG. 20. 

Since the retardation of tlie different DNA variants in the sieving matrix is a 
gradual process, separation of some homoduplexes and heteroduplexes may be 
achieved even outside the optimal denaturing conditions. The fom* lines in FIG. 21 
demonstrate the observed migration (scan nmnber) of two homoduplexes and two 

20 heteroduplexes in heterozygous lymphotoxin alpha samples (n = 1 1) analyzed at six 
different denaturing temperature conditions. The two lower lines represent the 
homoduplexes and the tsvo upper lines the heteroduplexes. At low denaturation, 
separation of the two heteroduplexes is achieved, as the heteroduplexes have 
thermodynamics with a lower stability than the homoduplexes. The standard 

25 deviation (e.g., variation in scan number) caused by the local differences in 

temperature in tlie capillary chamber is depicted in FIG. 21 . The result resembles 
the "S" shape observed in perpendicular DGGE, 
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The target lytnphotoxin alpha samples were then analyzed under different 
denaturing conditions. FIG. 22 A shows two electropherograms of the sarae sample, 
analyzed in the same run and in different capillaries under constant denaturant 
conditions. Since the capillaries are exposed to different temperatures (z. e., local 
5 temperature variation in the capillary chamber, see FIG. 20), the separation of 
homoduplexes and heteroduplexes differs in the two capillaries. Such capillary to 
capillary irreproducibility clearly prevents usage of a commercial capillary array 
system for CDCE type analysis. The PGR fragments shown in FIG. 22B were 
analyzed under cychc temperature conditions, from 48 °C to 52 °C in 15 cycles. The 

1 0 delay at each boxmdaiy temperature was set to 30 seconds. With the cycling 
denaturing conditions, identical separation patterns were repeatedly obtained in 
different capillaries. 

The reproducibiUty of cycling temperature capillary electrophoresis (CTCE) 
was evaluated by measuring the migration times of all peaks in a KRAS construct 

15 standard (see Examples 1 and 2). Separation of standard peaks and box plots, 

representing 649 samples standardized to migration time, is illustrated in FIG. 23. A 
normaUzed 

migration time was calculated by dividing peak migration time (i.e., the scan 
number) by the sum of all migration times of all pealcs in the same capillary. 

20 Discussion 

Capillary electrophoresis on multi-capiUary instruments, with the use of 
partial denaturing conditions, represents a recent approach in DNA variation 
analysis. Although this method is based on an established melting gel theory, 
optimization of denaturing conditions has generally been necessary by computer 
25 prediction and experimental verification for each individual target sequence. 
Described herein is the evaluation of tlie use of alternative temperature control 
strategies in order to achieve optimal denaturing conditions for analysis on a 
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standard commercial 96 capillary instrument. PCR fragments, previously evaluated 
by different melting techniques were used for analysis. 

The melting characteristics of ttie target sequences were first calculated by 
the Poland algorithm. The resulting theoretical thermodynamic behavior of the 
5 fragments was compared to experimentally observed migration of homoduplex 
forms of each fragment. An acceptable correlation was found between the 
calculated 50% melting probability values for the simulation with respect to 
migration differences between the two separated alleles. Consequently, melting 
algorithms developed for use in DGGE could directly be utilized on the capillary 
10 platform. However, at probability levels other than 50%, tlie calculated values were 
not in agreement with experimental results for the different homoduplex variants. 
The Poland web page has an option for estimating fragments containing mismatches 
(i.e., heteroduplexes). These tlieoretical simulations did not correspond well with 
the observed retention of heteroduplexes by cychng temperature capillary 
15 electrophoresis (data not shown). Apparently, theoretical calculation of 
themiodynamics of heteroduplexes in capillary melting techniques and its 
correlation with experimental data needs further investigation. 

By the use of WinMelt, we calculated the 50% melting probability of all 
oncogenic mutations in KRAS exon 1, and a linear correlation was found between 
20 relative electrophoresis retention of fragments and the relative thermodynamic 

stability. This correlation was used to determine the relative migration of a variant 
. fragment prior to the actual experiment. Based on the above thermodynamic 
calculation, the relative positions of mutant homoduplexes when related to wild-type 
could also be predicted. The relative position of a homoduplex in an 
25 electropherogram could then be used for direct identification of the variant from a 
set of possible sequences. This been demonstrated with KRAS exon 1. 

A critical parameter for denaturing CE with multiple capillaries is the local 
temperature fluctuations within the capillary chamber. As shown in FIG. 20, 
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differences across different capillaries were measured with both constant 
temperature and a single-sweep temperature gradient As a consequence of 
inhomogeneous temperature, loss of separation was observed in some c^illaries 
(FIG. 22). The main reason for these fluctuations was probably a non-equal air 
5 circulation inside the capillary chamber. As a consequence of the inaccurate 
temperatm:e control, an attempt to equiUbrate the differences by introducing 
temperature gradients and temperature cycling profiles was made. 

It was expected that in order to assure equivalent conditions for all separated 
samples, the gradient should be symmetrical, i.e., samples should spend equal time 

10 at each temperature. Fast cycling should therefore better serve separation 

reproducibility among parallels. During construction of different temperature 
cycling profiles it was observed that the maximum achievable temperature ramp rate 
was approximately 0.1 K/s. By ramping the temperature between maximum and 
minimum temperatures in the cycles and taking into accoxmt the cycling time, 

1 5 smooth sinus-hke temperature profiles were achieved. The maximum rate of cycling 
was then obtained fi"om the temperature range and the temperature ramp rate. It was 
observed that too high cycling of temperature resulted in lower temperature 
amplitudes in the sinus curve. A maximum of 20 cycles spanning 3 K in 20 min 
yielded reproducible temperature curves. Secondly, another notable advantage of 

20 fast cycling is the elimination of problems due to changes in the optical response of 
the instrument. Most commercial instruments include precise optical elements 
(filters, lenses, pinholes, mounts, holders, etc.) that are susceptible to temperature 
changes. The usual impact of changing temperature on optics includes loss of 
aUgnment and focus, resulting in reduction of signal-to-noise. With fast temperature 

25 cycling, most optical elements, which are made of solid materials, do not follow the 
rapid changes. As a result, their temperature is equilibrated at an average constant 
temperature and maintained during the entire run. The reproducible performance of 
the instrument is therefore maintained during the temperature gradient operation. 



wo 2005/047541 



PCT/IB2003/005304 



-57- 

Thirdly, use of a cycling gradient further allows significant increase in sample 
tbroughput by application of multiple-injection workflow (Example 1). Most 
commercial capillary array instruments allow controlled temperature during the run, 
thus this approach of DNA variation analysis should not be instrument dependent. 
5 The use of denaturant gradients with Joule heating and temperature zone 

have been appUed to single capillary instruments (Gelfi, C. et al, 1996, 
BioTechniques, 21:926-932; Zhu, L. et al., 2001, Electrophoresis, 22:3683-3687; 
Schell, J. et al., 1999, Electrophoresis, 20:2864-2869). However, these strategies 
did not improve the sensitivity and specificity of the mutation analysis as compared 
10 to CDCE (Khrapko, K. et al.. Nucleic Acids Res., 22:364-369; Gelfi, C. et al., 1997, 
Electi-ophoresis, 18:724-731). hi the present setup, the use of gradients and cycUng 
in the multi-capillary system was adapted for an entirely different reason: The 
original purpose of the gradient was to ehminate local temperature differences in the 
capillary chamber. Unexpectedly, improved separation in each capillary was 
15 observed (FIGS. 22 and 23). However, the benefit of denaturing gradients applied to 
slab gel techniques (TGGE, TTGE, and DOGE), with regard to robustness, appeared 
as a secondary effect of a cycled temperature gradient. In both gradient slab gels and 
CTCE, fine-tuning of denatiu-ant conditions is not needed, as is the case in CDGE 
and CDCE. Sensitivity and specificity reported for CDCE is preserved in the CTCE 
20 setup. This is because the separation time in CTCE (i.e., the time each capillary 
spends at the optimal denaturing temperature) is comparably much longer than for 
single gradients. 

In a common configuration of DGGE, analysis is performed on a slab gel 
with defined gradient of chemical denaturant perpendicular to the electrophoretic 
25 direction. The resulting "s-shaped" diagram is used to locate optimal denaturing 
conditions. In a capillary format, a similar analysis can be performed by comparing 
separation profiles of samples separated at different denaturing conditions. The 
results presented in FIG. 21 provide information on the temperature interval in 
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which separation of homoduplexes and heteroduplexes is achieved. Furthermore, 
this can be used in settings were different target fragments are to be analyzed in the 
same run. By adopting a compromise with regard to the denaturing conditions for 
individual fragments, samples differing by up to 6 K in thermodynamics of the low 
5 melting domain could be analyzed in parallel capillaries (data not shown). 

Additionally, running the fragment under multiple different denaturing conditions 
helps to discriminate artifact peaks formed during PGR or spikes in electrophoresis. 
Their migration will not be significantly altered by small changes in denaturing 
conditions, and can thus be excluded as aberrant (Bjorheim, J. et al, 2002, J. Sep. 

10 ScU 25:637-647). 

To demonstrate reproducibiUty of the temperature cycling, 7 subsequent runs 
of an artificial KRAS construct exhibiting 5 peaks (representing 5 different KRAS 
exon 1 variants) were performed. The migration times of all peaks were normalized 
to compensate for offset in individual traces (e.g., due to differences in 

15 electrophoretic current). Variance of the normalized migration times is shown ia 
FIG. 23. The reproducibility of the method is demonstrated by the mean and 
standai'd deviation iq the box plot transposed over a representative electropherogram 
of the KRAS constract. It can be clearly seen from this figure that the individual 
boxes are well-spaced. This indicates a complete baseline resolution of all peaks at 

20 any given combination of runs or capillaries. 

Denaturant capillary electrophoresis has previously been demonstrated to 
have high sensitivity and high throughput potential in DNA variation analysis. With 
CDCE, time-consuming optimization has to be expected as well as lower 
reproducibility compared to cycling of temperature in denaturant CE. Consequently, 

25 mutations in target sequences have been detected, but their identity had to be 

revealed by other means (Guldberg, P. and Guttler, F., 1993, Nucleic Acids Res., 
21:2261-2262). Analysis with cycling temperature exhibits a high robustness and 
reproducibility allowing direct identification of mutations based on peak migration 
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times in electropherograms. The results presented in this study can be directly 
transferred to any multi-capillary sequencing instrument with air-ventilated chamber 
without any need for hardware changes. In theory, all fragments previously analyzed 
with denaturing 

slab gel techniques can be analyzed with cycling temperature in a denaturant 
capillary electrophoresis setup without any need for complex PGR procedures. In 
conclusion, the theory of well-known denaturant slab gel techniques can be 
transferred to the multi-capillary format applied to denaturing CE. This opens a 
possibility for fully automated analysis of DNA variation in large sample sets for 
high throughput screening. 

Example 6, Simultaneous allele separation of DNA vaHants in twelve different 
fragments, with large differences in thermodynamics, by high throughput cycling 
temperature capillary electrophoresis. 

Described herein is an analytical tool for parallel defection and screening of 
15 nucleotide variation (mutations and single nucleotide polymorphisms) in DNA 

fragments. Twelve different alleles with known single base DNA variation sited in 
eleven PGR fragments with theoretical melting temperature ranging from 69 to 78 
*^G were amplified in 96-well plates by a conmion protocol. Following PGR, 
separations of alleles were achieved by cycling temperature capillary electrophoresis 
20 (GTGE) in a standard DNA sequencing instrument. By applying denaturiag 

temperature in cycles, covering the theoretical melting temperatures of all fragments, 
separation between variants and canonical alleles were achieved in all samples in 
parallel runs. Smaller temperature cycles resulted in loss of separation between 
homoduplexes in fragments with theoretical melting temperatures not included in 
25 the cycling temperature gradient. On the other hand, separation between variants 
was greater with more focused cycling gradients than with a general large 
temperature cycle gradient. By adopting a compromise with regard to the denaturing 
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conditions for individual fragments, samples with large difference in 
thermodynamics were analyzed in parallel capillaries with CTCE. The method 
opens a possibiUty for one step detection and screening of DNA variants in 96 
different fragments without need for individual analysis design and optimization. 
5 By introducing large temperature gradients, separation of DNA variants with 

large differences in thermodynamics is achieved by CTCE. This allows analysis of 
different DNA fragments in the same run. The analysis is performed in a 96 
capillary instrument with four channels, allowing simultaneous analysis of up to 384 
different samples within 40 minutes. Further, the fragments can be scanned.for 
10 unknown DNA variation as previously reported for melting gel techniques. For 
demonstration purposes, twelve DNA variants with 9 K difference in 
thermodynamics were analyzed imder various denaturing conditions in an automated 
and high throughput system with the aim of separation all variants in parallel runs. 

DNA samples and PGR 

1 5 Genomic DNA was extracted from anonymous blood donor samples by 

standard DNA extraction kit (Qiagen Inc., Valencia, CA, USA). Primers were 
designed with PrimerS (Whitehead Institute, Cambridge, MA) software for the target 
sequences, using standardized parameters, thus all fragments could be amphfied 
with the same PGR conditions. Twelve different samples were analyzed in 

20 replicates of eight (12 DNA fragments x 8 wells = 96 wells) in 96-well plates. 

AmpUfication was performed on a PTC-200 thermal cycler (MJ researcli, Waltham, 
MA, USA), by mixing 50 ng genomic DNA with 25 \iM of each dNTP (Perkin 
Ehner, Oslo, Norway), lOX Tag buffer, 1 imit Tag and 5 pmol of each primer 
(MedProbe, Oslo, Norway) in a final volume of 25 (iL. The cycUng parameters were 

25 35 cycles of denaturation for 30 seconds at 94 °C, annealing at 56 **C for 30 seconds 
and elongation at 72 for 1 minute, followed by elongation step of 72 "^C for 10 
minutes. 
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Target sequence identification 

Fragments were design for analysis of DNA variation identified by NCBI 
SNP reference number rs, 1) Breast Cancer, Type 1 (BRCA 1) rs799923, 2) BRCA 1 
rsl6940, 3) 5, 10-Methylenetetrahydrofolate Reductase (MTHFR) rslSOllSS, 4) 
5 Opsin cone pigments (OPNl) (base variation between long and medium wave genes 
at base pair 78173 in Genbank sequence AC092402), 5) MTHFR rslSOl 131, 6) 
MTHFR rs2274976, 7) Cystathionine-Beta-Synthase (CBS) rs234706, 8) Nad(P)H 
Dehydrogenase, Quinone 1 (NQOl) rsl 800566, 9) Dihydropyrimidine 
Dehydrogenase (DPYD) rs3918290, 10) DPYD (-10 base pair firom rs3918290), 11) 
10 DPYD rslS01265, 12) Cytotoxic T Lymphocyte-Associated 4 (CTLA-4) rs5742909. 

Theoretical thermodynamics 

Theoretical melting profiles of the different target fi-agments were obtained 
by applying Poland algorithm OPoland, b., 1974, Biopolymers, 13:1859-1871) 
throu^ the iise of WinMelt computer program (MedProbe, Oslo, Norway). The 
15 optimal separation temperature proposed by these programs was adjusted on the 
basis of the urea concentration in the matrix. For each molar increment of urea, the 
temperature was lowered approximately 3.2 K (Fodde, R. and Losekoot, M., 1994, 
Hum. Mutat., 3:83-94). 

CE 

20 PCR samples were analyzed in a standard 96-capillary DNA analyzer 

(MegaBACE™ looo DNA Analysis System, Amersham Bioscience, Sunnyvale, 
CA, USA). The distance jfrom the anode to the detector was 40 cm. Linear 
polyaciylamide (MegaBACE LP A) containing 7 M urea was replaced in capillaries 
prior to each run. Standard coated capillaries were used. PCR products were loaded 

25 into the capillaries by electro kinetic injection at 200 V/cm for 12 seconds. The 
electrophoresis was carried out at a constant field of 200 V/cm. Laser-induced 
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fluorescence was used with excitation at 488 lun (blue laser) and detection of 
emission at 520 nm (FAM channel). 

Denaturing conditions and temperature control 

In each run, a 96-well plate with the same distribution of fragments was 
5 analyzed xmder various denaturing conditions. The denaturing temperature in the 
capillary chamber was either cycled with large amplitudes covering the theoretical 
thennodynamics of all selected fragments (13 K) or with smaller amplitudes down to 
3 K. After the cycling, the temperature in the chamber was reduced to 40 "^C. The 
denaturing conditions, e.g., the temperature profiles, were controlled through the 

10 Instrument Control Manager (ICM) software package (Amersham Bioscience). The 
MegaBACE™ Sequence Analyzer View and Edit software program (Amersham 
Bioscience) was used to measurer the migration times of all peaks. Three cahbrated 
temperature sensors were placed in close proximity to the capillaries at the same 
vertical level, but at different horizontal locations. One temperature sensor was 

15 placed in tlie front part, one in the middle part and one at the back of the capillaries. 
The temperature recording from the three thermocouple probes (K type, Sivert Max, 
Oslo, Norway) was recorded using an analogue computer board (lotech, Inc., 
Cleveland, Ohio, US A) at a rate of 1 Hz. 

Results and discussion 

20 Theoretical melting profile analyses of the selected target fragments with 

attached GC-clamp were first performed. The thermodynamics of the canonical 
DNA sequences of all twelve fragments are illustrated in FIG. 24. The 
themiodynamics of the low melting domain of tlie CTLA-4 sequence showed a 
theoretical melting temperature, 68.7 ^C. The selected sequence from the CBS gene 

25 had the highest thermodynamics of the low melting domain, 78.1 ^C. All other 
fragments had theoretical melting temperatures between these extremes. Variant 
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sequences were also analyzed and demonstrated a theoretical difference in 
thermodynamics as compared to wild-type sequences (data not shown). The 
selected fragments thermodynamics are a representative selection of ampUcons 
designed for melting gel analysis. On average, fragments designed on random DNA 
5 will display thermodynaihics between 70-75 "^C. This is exemplified by computing 
the thermodynamics of a random 20000 base pair in chromosome 13. The average 
thermodynamic was found to be 68.9 **C per base pair, with iso-melting domains in ■ 
the range of 100 to 200 base pairs. Thus, theoretically, 96 different fragments can be 
constructed with thermodynamics similar to that of the fragments used in this study. 

10 After evaluation of the thermodynamics, creation of a 96-well plate with 

replicates (n = 8), of the different PGR amplified fragments were electrophoresed 
under various denaturing conditions. When temperature in the capillary chamber 
was cycled between 59-47 ""C (a denaturing temperature interval that covers the 
thermodynamics of all of the fragments) five times, separation of all twelve 

15 fragments were achieved repeatedly in parallel runs and in different capillaries (FIG. 
25). Separation of homoduplexes and heteroduplexes was performed within 40 
minutes for all samples. As seen by the base line separation between the 
homoduplexes (peak 1 and 2) variant homozygous samples can directly scored in 
electropherograms by use of an internal standard, with a different label {e.g., 

20 TAMRA). This is in clear contrast to methods like heteroduplex analysis (HA) and 
dHPLC,, where separation of the homoduplex species has not been reported, except 
for large deletion/insertions (Cremonesi, L. et al, 2003, Br. J. Haematol, 
121 :173-179; Kozlowski, P. and Krzyzosiak, W., 2001, Nucleic Acids Res,, 29:E71). 
The migration times of all peaks of all fragments in a 96-well plate analyzed with 

25 the above-described denaturing conditions were measured. The values found for the 
second homoduplex and heteroduplexes were normalized to the migration of the 
first homoduplex peak. In FIG. 26 the mean migration times of peak maximimi for 
the homoduplex and heteroduplex peaks is shown. The brackets indicate one 
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Standard error of the mean in FIG. 26. Baseline separation of all homoduplexes was 
achieved under these cycling conditions. Occasionally, some of the heteroduplex 
peaks co-migrated with the large temperature gradient, due to complete melting of 
the low melting domain of both heteroduplexes. Incorporated in the figure is an 
5 electropherogram presenting analysis of firagment number 6. This is the first time 
that separation of twelve different alleles has been achieved to such an extent that 
direct genotyping is possible. Thus, comparing this method with exciting DNA 
variation discovery methods, dHPLC, HA or SSCP renders this adaptation of DGGE 
on to a regular capillary DNA sequencing instrument the most powerful tool in DNA 
10 discovery. 

Constant denaturing capillary electrophoresis results were dependent on very 
precise temperature control to be able to reproducibly separate alleles 
(Li-Sucholeiki, X. et al, 1999, Electrophoresis, 20:1224-1232). With this in mind, 
measurement of the denaturing temperature at dijBferent locations within tlie capiUary 
15 chamber was performed during electrophoresis. FIG. 27 shows the observed 

denaturing temperature in the capillary chamber with a preset cycling temperature 
gradient ranging from 59-47 in 5 cycles. Due to the large difference between 
minimum and maximum temperature and taking into account the cycling speed, 
temperature within tlie capillary chamber did not reach the extreme preset 
20 temperatures during electrophoresis. Further the temperature at different locations 
in the capillary chamber was incongruent, especially at high denaturing 
temperatures. However, this discrepancy in microenvironment between capillaries 
with in the chamber did not have any significant effect on the separation (FIG. 26). 
The observed difference betv^^een maximum and minimum denaturing temperature 
25 within the capillary chamber was 8 K. As previously reported, fast cycling of the 
temperature may lead to narrowing of the temperature amplitude (Example 5), with 
loss of separation in some capillaries as a result. By applying large temperature 
gradients or temperature cycles during electrophoresis, the firagments are exposed to 
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sub-optimal denaturing conditions for a longer time as compared to more focused 
denaturing conditions. Nevertheless, as long as the maximiim denaturing 
temperature is below the denaturing temperature.of the GC-clamp, separation of 
variants is preserved during electrophoresis (data not shoAvn). By analyzing 
5 ' fragments at increasing temperature gradients, or increasing the average gradient 
temperature, specific "melting" characteristic of the low melting domain can be 
observed. Raking the peaks and subtracting the elution time of each species gives 
six different melting curves. FIG. 28 demonstrates the relative distance between 
peaks in fragment nmnber 1 (BRCAl) at different denaturing temperatures. The 
1 0 optimal denaturing temperature of this fragment was calculated to be 48.6 ^'C, a 
value that correspond well with separation maximum observed in FIG. 28. 
Separation of fragments at other denaturing temperature can also be expected as seen 
in the figure, and this quality allows compromising of the denaturing conditions in 
parallel analysis of fragments with different thermodynamics. Noteworthy is the 
1 5 melting of the heteroduplexes, which are, shifted a couple of degree lower, than 
homoduplexes, due to the miss-match in double stranded DNA (FIG. 28). 

Cycling denaturing temperatures with smaller amplitudes was also 
investigated. Experiments demonstrated that with smaller denaturing cycles only 
fragments with theoretical melting temperature corresponding to the denaturing 
20 conditions separated during electrophoresis (data not shown). On the other hand, the 
separation between homoduplexes and heteroduplexes were better than with large 
cycling gradients. FIGS. 29 and 30 illustrates the observed denaturing temperature 
in tiie capillary chamber with a programmed gradient from 47 to 55 "^C, and the 
corresponding separation of fragments during electrophoresis, respectively. A large 
25 discrepancy between progranmied and observed temperature was observed, when tlie 
start temperature was at the low extreme. This can be explained by the heat capacity 
of the instrument. With this setup the homoduplex in fragment 7 did not separate 
out, due to low maximum denaturing temperature. Nevertheless, separation of 
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heteroduplexes was observed in all runs. The observed temperature during 
electrophoresis and the measured separation of alleles correspond well with the 
thermodynamics results described in FIG. 24. 

Several analytical tools are designed for diagnostic analysis of DNA 
5 variation. One increasingly popular approach is dHPLC, which is based on similar 
separation principles, as the method reported here, were melting of double stranded 
DNA is separated in colvmms. By comparing the results reported herein and the 
Uterature one has to draw the conclusion that CTCE is superior to dHPLC by several 
means. The argument is described in the following paragraphs. 
1 0 First, the cost of a dHPLC system is equivalent to a 48 capillary DNA 

sequencing system (MegaBACE 500 TM) were the obviously fact is that the CE 
instrument can analyze 48 samples at the time compared to one in the dHPLC 
system, when using only one channel of the DNA sequencing instrument. 
Acknowledging that DNA sequencing instruments in general have four channels, the 
15 sample throughput becomes hundred fold, that of the "high throughput" dHPLC 

(Xiao, W. and Oefiier, P., 2001, Hum. Mutat, 17:439-474). Furthermore, previously 
reported sample throughput of CTCE is limited to approximately 12000 samples per 
24 hours in a 96 capillary instrument (Example 1). This is in gross contrast to the 
alleged high throughput of dHPLC, with a capacity of 144 samples m 24 hours. 
20 Second, by implementing this assay one a regular DNA sequencing 

instrument, three application can be run on one instrument (DNA sequencing, 
fragment analysis and DNA variations assay) while only one assay, detection of 
DNA variation, can be run on dHPLC. hi addition, given that DNA variations are 
detected by dHPLC a need for a second instrumentation {i.e., DNA sequencing 
25 instrument), or other means are needed to determine the type and position of the 
DNA variant (Hoogendoom, B. et al., 1999, Hum. Genet, 104:89-93). 

Third, the sensitivity defined as ratio of mutant versus wild type is reported 
to be in the 5 to 10% range for dHPLC (Wolford, J. et al, 2000, Hum.GeneL, 
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107:483-487), while denaturant capillary electrophoresis has been reported down to 
• 0.1% (Li-Sucholeiki, X. and ThiUy, W., 2000, Nucleic Acids Res., 28:E44; Khrapko, 
K. et al, 2001, Methods Mol Biol, 163:57-72). Consequently when analyzing 
heterogeneously sample like tumor tissue or pooled DNA, low mutant fraction 
5 cannot be detected by dHPLC. 

Fourth, as demonstrated herein implementing firm knowledge of DNA 
thermodynamics and applying denaturing cycling gradient thereafter, any variation 
witliin the fragment can be sq)arated on heteroduplex level and a large fraction on 
the homoduplex level. While dHPLC has to fine-tune the denaturing temperature 
10 for each variation to be analyzed, whether it is in the same fragment or in a different 
amplicon. 

Last, there are no reports on separation of homozygote variants by dHPLC; 
save for large deletion on one allele, hence homozygote mutation or polymorphisms 
will not be detected. On the contrary, CTCE have clearly demonstrated baseline 

1 5 separation between homoduplexes and heteroduplexes in the results shown herein 
and in other reports (Example 2; Bjorheim, J. et al, 2003, Mutat Res., 526:75-83; 
Kristensen, A. et al, 2002, BioTechniques, 33:650-653). Clearly the arguments 
above render CTCE a superior with regard to cost, sample throughput, specificity, 
sensitively and simplicity (one gradient fits all) as compared to dHPLC. 

20 The analytical tool reported herein allows parallel analysis of different PCR 

products with up to 9 K difference in theoretical thermodynamics by CTCE. 
Ninety-six fragments are analyzed within 40 minutes, and sequence variants can be 
directly identified when control samples are available. If PCR primers are designed 
with almost identical T^, different DNA sequences can be amplified simultaneously 

25 in 96-well plates, and the plates can thereafter be directly transferred to CTCE 

analysis. Since the fragments are exposed to cycling denaturing temperature during 
electrophoresis, new series (e.^., new 96-well plates) of fragments can be repeatedly 
be injected in the capillaries colimms as previously reported (Example 1). The 
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sample tooughput of the method in diagnostic analysis is therefore, for any practical p 
reasons, beyond obtainable number of patient samples in most settings. Further, 
different DNA sequences can be included in designed diagnostic 96-plates, relevant 
for dijfferent clinical settings and hypotheses, making the method amenable for a 
5 huge research community. Finally, the approach can easily be adapted to a fiilly 
automated set-up with the Caddy™ - robotic sample plate manipulator used in 
conjunction with the MegaBACETM, which allows continuous analysis of samples 
without need of an operator. 

Example 7. Analysis of DNA variation by denaturaiit capillaiy electi-ophoresis, 
10 For all of the applications described below, alleles (DNA variants) are 

amplified in the same polymerase chain reaction (PGR) reaction by the same primer 
pair, subsequently followed by allele separation by electrophoresis. 

PGR 

To efficiently analyze DNA, the number of copies in a sample first was 
15 mcreased. This is performed through the well-characterized polymerase chain 

reaction. PGR does not only increase number of copies of DNA exponentially, but 
also restricts the amplification to specific target sequences determined by the 
primers. Specific PGR protocols are common in the art. 

Electrophoresis 

20 Gapillary Electrophoresis (GE) was first described in 1967. During GE, a 

sample is separated into its components as it migrates through a capillary under the 
driving force of an electric field. Separation is typically performed in a long (1 0 to 
100 cm), narrow (10 to 100 |xm), electi-olyte-filled, fiised siHca capillary, hi the GE 
systems used to study DNA variants, a sieving matrix, e,g., chains of polymer, are 

25 used as separation media. 
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DNA MELTING GEL ANALYSIS 
Theory Of DNA Thennodyaamics 

In 1974, Poland proposed an algorithm able to calculate the melting 
probability of thousands of nucleotides in dsDNA (see above). Based on nearest 
5 neigihbor correlation in q)ecific sequence macromolecules and basic principles of 
DNA thermodynamics, the algorithm stated that dsDNA melts to single stranded 
DNA (ssDNA) when exposed to sufficiently high temperatures and/or chemical 
denaturants (e,g,, formamide and m^ea). Specifically, the length of the DNA 
fi-agment and the nucleotide sequence within the firagment defines the melting 

1 0 temperature at which each base pair (bp) of a DNA duplex is in perfect equiUbrium 
between the denatured'and helical state. Importantly, DNA variants differing by 
only one base will reveal different melting properties based on the sequence 
variation. Differences in thermodynamics between differing DNA sequence such as 
base substitutions are exploited in melting gel techniques. 

15 The need for approximations may rely on the fact that an exact algorithm 

required computer time proportional to N2, were N equals the number of bp in the 
target dsDNA. Computer software programs such as SQHTX, Melt87, WinMelt and 
the Poland internet web site calculate the melting profile of defined homozygous 
dsDNA fragments based on different related algorithms. Typically, stretches of 

20 double-stranded DNA consist of iso-melting domains with local thermodynamics 
covering a few hundred bases. The part of a firagment that melts at low denaturant 
exposure is known as the low melting domain, whereas the part that melts at high 
denaturant exposure is called the high melting domain (FIG. 31). The terai 
"melting" refers to the change in the structure of DNA from an orderly helix to a 

25 disordered structure without base paring. The programs calculate the midpoint 

temperature at which each bp is at 50/50 equilibrium between the helical and melted 
states. Since the target sequence must be in the low melting domain, not all 
fragments will reveal a melting profile directly suitable for melting gel analysis. If 
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the target sequence is part of a GC-rich area, it will often be part of a high melting 
domain. The goal of melting profile analysis of DNA fragments with computer 
programs is to select and manipulate the target sequence so that the region of interest 
is in the low melting domain. This often requires the proper placement of an 
5 artificial high melting domain. For a detailed description of firagment construction 
and GC-clamp attachment see (Example 5; Bjerheim, J. et aL, 2003, Mutat. Res.^ 
526:75-83). In theory, almost any DNA sequence can be analyzed with melting gel 
techniques as long as the base change results in sufficient differences in 
themiodynainics between alleles and that the target sequence is located in a low 
10 melting temperature domain adjacent to a high melting temperature domain (FIG. 
31). ~ 

Denaturant Gradient Gel Electrophoresis and related techniques 

Different melting gel techniques are named after the combination of 
denaturants in the gel (slab-gels) or matrix (capillary platform) and the platform on 

15 which the analysis is performed. Denaturant gradient gel electrophoresis (DGGE) 
has been tlie most frequently used meltiag gel techniques and was first described by 
Fisher and Lennan. DGGE has successfully been applied to detect mutations in 
genes including PKU, TP53, KRAS, NRAS, HRAS as well as genes involved in 
cystic fibrosis (Fischer, S. and Lerman, L., 1980, Proc. Natl Acad, ScL USA, 

20 77:4420-4424; Eileen, H. et al, 1996, Hum. Mutat, 8:19; Nedergaard, T. et al, 
1997, hit, J. Cancer, 71:364; Wang, J. et al, 2000, Mol. Genet Metab., 70:316). 
All other melting techniques are variants of this method. During Constant 
Denaturant Gel Electrophoresis (CDGE) the fragments are subjected to constant 
denaturing conditions in order to achieve better separation between alleles. CDGE 

25 has been used to detect base pair changes in defined target sequences such as TP53, 
BRCAl, RBI, hMSH2, MSH6, RET and HPRT (Smith-Scrensen, B. et al, 1992, 
Mutat Res,, 269:41; B0rresen, A. et al, 1991, Proc, Natl Acad. Set USA, 88:8405; 
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Hovig, E. etal, 1992, Genes Chromosomes. Cancei\ 5:97; Bfioresen, A. et al, 1995, 
Hum. Mol. Genet,, 4:2065). On the other hand, the use of CDGE regulars carefixl 
casting of denaturing gels and control of denaturing temperature. If the denaturing 
conditions are slightly off scale separation will be lost. Other modifications of 
5 DGGE have focused on changing the temperature surromiding the gel and have 
resulted in methods including Temporal Temperature Gel Electrophoresis (TTGE; 
Yoshino, K. et al, 1991, Nucleic Acids Res., 19:3153) and Temperature Gradient 
Gel Electrophoresis (TGGE; Riesner, D. et al, 1989, Electj'ophoresis, 10:377). The 
denaturing conditions are therefore gradually changed during electrophoresis. The 

10 two main advantages of TTGE and TGGE compared with DGGE and CDGE are: 1) 
The simpler acrylamide gel casting procedure and 2) the potential to analyze 
fragments containing more than one isomelting domain. Both TTGE and TGGE 
have been employed in the analysis of KRAS, PTEN, mitochondrial DNA, TP53 
and HLA. All of the above mentioned slab-gel methods have to some extent been 

1 5 extended to the capillary foraiat; 

DENATURANT CAPILLARY ELECTROPHORESIS 

The modification of DGGE to the capillary platform has reduced the 
separation time of mutants, increased s^aration efficiency, and enabled the 
detection of low-frequency mutations in human cells and tissues with mutant 

20 fraction down to 10'*^. 

The application of Constant Denaturant Capillary Electrophoresis (CDCE) 
was developed on a laboratory assembled instrument and later adapted to 
commercial available capillary DNA sequencing instruments. In CDCE the optimal 
denaturing conditions for separation of DNA variants are deteraiined by fine-tuning 

25 the temperature by re-analyzing the same sample at different temperatures. Within a 
certain temperature range, the position of the melting equilibrium and thus the 
average electrophoresis mobility of each mutant is different, which allows sequences 
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containing single base pair point mutations to be separated. The CDCE method has 
been appUed to analysis of DNA variation of several fragments including 
mitochondrial DNA APC, HRAS, KRAS and TP53. A further improvement of 
CDCE by modifying the instrumentation to allow two-point detection and automated 

5 fraction collection has demonstrated detection of mutations in samples with mutant 
fractions as low as 10"^. Despite the sensitivity and specificity of the method, few 
laboratories have applied this technique. This may rest upon the fact that polymer 
replacement and sample injection is performed manually, so the analysis is labor 
intensive and dependent on an instrument operator. Conversion of CDCE onto 

0 regular capillary DNA sequencing instruments has resulted in automation and 

standardized protocols. Compared to conventional CDCE in which each sample has 
to be loaded manually, automated CDCE analysis can be performed in a single or 

- multi capillary instrument. The automation of CDCE allows for rapid analysis of a 
large number of samples over a short period of time and virtually any capillary 

5 sequencing instrument with a temperature controller can be used. The method was 
first described with use of an ABI 310 Genetic Analyser, in which up to 48 samples 
can be analyzed by this method without any need for operator intervention. The 
method has also been adapted to a regular 96-capillary DNA sequencing instrument, 
the MegaBACE™i000, allowing for analysis of 96 within 40 minutes. 

3 The introduction of a temporal temperature gradient within the capillary was 

produced by a voltage ramp during electrophoresis. Li this technique the 
temperature increase within the capillary due to the Joule heating during 
electrophoresis results in separation of alleles (Gelfi, C. et al, 1996, BioTechniques, 
21 :926-932). A different approach was later taken by modifying the CE system with 

5 a thermostatic liquid surrounding the capillary to allow for computer generated 
temporal temperature profiles between 25 ""C and 70 ""C. A temperature gradient 
was also used in an automated 96-capillary DNA sequencing instrument to 
discriminate between mutations. Samples were analyzed with a temperature 
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gradient, starting above and declining beneath the optimal separation temperature, 
by control of the instrument software. This technique proved to be robust in that the 
gradient compensated for the temperature differences within the capillary chamber, 
moreover each capillary passed through the optimal separation temperature around 
5 the theoretical melting temperature for the analyzed fragment. There was no need 
for time-consuming optimization procedures once the gradient has been established 
around the theoretical melting temperature of the fragment. Further improvement of 
the gradient concept was accomplished by cycling the gradient several times around 
the optimal separation temperature. This is based on the observation that during the 
1 0 cycling gradient, each capillary will pass through its optimal melting condition for 
the fragment of interest several times. This method is more robust than the CDCE 
when applied to multi-capillary instruments due to temperature differences within 
the capillary chamber. 

Mutation analysis 

15 Detection of genetic changes is of central importance in molecular genetics 

and cancer research since alteration of just a single base pair in a DNA sequence can 
lead to change in cellular behavior. Knowledge of genotype markers, related to 
tumor behavior and patient prognosis, may lead to more effective preventive 
measures, prediction of prognosis and better personal medical treatment. 

20 Furthermore, the possibility of detecting mutations m a low fraction can be 

important for early detection of malignant diseases, detection of remaining cells after 
surgery and possibly for prognosis and outcome of the disease. 

In two reports, denaturant capillary electrophoresis (DCE) was applied to the 
detection of mutations in KRAS exon 1 and exons 5-8 of TP53 gene by use of an 

25 ABI 310 Genetic Analyzer, with a sensitivity (mutant/wild-type) down to 1 0"^ 
(Bj0rheim, J. et al^ 2001, Tumour, Biol, 22:323-327; Bjorheun, J. et al, 1998, 
MutatRes.^ 403:103). The same application was also demonstrated on a modified 
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96-capillary DNA sequencing instrument, the MegaB ACE™i OOO. The need for 
modification of the instrument rested upon the need for accurate temperature control 
in all 96 capillaries, so an extemal heater was fitted to the capillary array. To 
circumvent the use of an extemal heater a temperature gradient controlled by the 
5 instrument software was later applied to the sequencer. 

This modified technique was used to successfully separate mutants in exon 8 
of the TP53 gene in all 96 capillaries (Kristensen, A. et al, 2001, BioTechniques, 
33 :650~654). In a reconstruction experiment the sensitivity to detect mutations 
alleles in a wild type background was 0.4 %. To further improve the reproducibility 

10 of the separation in a multi capillary instrument, the element of cycling temperature 
during electrophoresis was introduced. By changing the temperature aroimd the 
optimal separation temperature, almost equal denaturing conditions in all capillaries 
and between runs was obtained. As a result, highly reproducible denaturing 
conditions were observed, which resulted in detection and identification of all 

15 twelve different KRAS exon 1 mutations (Example 2), Furthermore, we have 
applied DCE (CDCE, TGCE and CTCE) analysis of mutations in several genes 
(KRAS, TP53, NRAS, PTEN and BRAF). The figure below demonstrates detection 
of mutants in samples with different mutant firaction, by DCE (FIGr 32); 



Single Nucleotide Polymorphisntis (SNP) 

20 A SNP is defined as base substitution with a frequency above 1% in the 

population of interest. When two randomly chosen alleles are compared, sequence 
variations are expected to be found every 1000 bp. Consequently, one would expect 
many sequence differences when analyzing larger populations as compared to the 
canordcal sequence. In addition, SNPs are typically bi-allelic and are therefore 

25 becoming popular to discriminate between alleles or haplotypes in medical and 
population genetic studies. Some of these sequence variants may affect functions 
such as immune response, metabolism, DNA repair, or control of cell death and 
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n 

division. Finally, some variants may be regaided as SNP in one population aad as a 
mutation in another cohort. 

An approach utilizing DCE for SNP analysis is described herein. Target 
sequences were PGR amplified followed by allele separation by capillary 
5 electrophoresis. The genotypes of the individuals SNPs were scored based on 
co-migration to an internal standard. The firagment specific intemal standard was 
made from a sample with both variants (alleles) of the SNP and was labeled with a 
different fluorophore than the sample (FIG. 33). 

Due to the presence of the intemal standard in all the electrophoretic nms, 
10 the specificity of the SNP analysis is retained at an excellent level because allele 
separation of the standard is required before genotypes can be determined. 

Microhaplotypes 

A Microhaplotypes is in general defined as the presence of several SNPs 
within a short fi-agment of DNA (100-1000 bp). The theoretical number of possible 

1 5 alleles is defined as 2n, were n is number of SNPs in the sequence analyzed. 

Possible combinations of genotypes are given by the sxmi of nimibers fi-om zero up 
to numbers of alleles (e.g., 4 = alleles, genotypes = 4+3+2+1). Few methods are 
able to analyze microhaplotypes, and the protocols are in general labor intensive and 
not easily automated. Consequently, DCE has been optimized for microhaplotyping 

20 on some jfragments by use of standard DCE protocol, e.g. , by PCR and 

electrophoresis. DCE w£is used to identify multiple SNPs within a defined target 
sequence based on theoretical thermodynamics and migration of variant fi-agments 
(Krisetensen, A. et al^ in press). This target sequence is located in the 
Ataxia-telangiectasia (ATM) gene, in which three SNPs are reported. 

25 Fluorophore-labeled PCR products were analyzed without any post PCR steps on a 
standard 96 capillary- sequencing instrument MegaBACE™1000, 3677 samples 
were analyzed for microhaplotypes, to validate the method, and further to estabhsh 
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microhaplotype frequencies in the target populations. This direct SNP detection and 
haplotype determination based on co-elution have added one more level of 
information to the SNP analysis. In an ongoing study of familial breast cancer, 
microhaplotying is used to determine haplotype sharing in BRCA 1 and 2 genes and 
5 has efficiently determined microhaplotypes in several thousand individuals. 
Exemplified microhaplotypes from two adjacent SNPs located on chromosome 
1 8q2 L 1 are shown in FIG. 34. Alleles were separated by DCE in a standard DNA 
capillary sequencing instrument and genotyped according to an internal standard 
having allele 2 and 4 (FIG. 34). Electropherograms of selected samples are presented 
10 in FIG. 35. 

With only one PGR reaction and one electrophoresis two neighboring SNPs 
can be determined with the additional information of microhaplotype. The 
application has therefore been improved from labor intensive and low capacity to a 
straightforward and high-throughput method for microhaplotype detemiination. 

1 5 Gene Copy Number 

Pooling of samples is an obvious approach when analyzing SNPs in large 
populations. Described above is a method where the population frequency of SNP 
can be efficiently detected, and their allele frequencies accurately measured 
(Example 3). In this example, 8000 alleles were analyzed for selected SNPs. There 

20 was no need for any correction of the signal of the separated alleles because both 
alleles are labeled with the same fluorophore and because fragments of the same 
length with the differences of one base pair will pass the detector with the same 
velocity. Thus, quantum yield will reflect DNA copies entered in the PGR. 
Additionally, no preferential ampUfication of alleles could be observed when areas 

25 under the peaks were measured in heterozygous individual samples. The ratio of 
alleles was essentially 1:1. The result of the study gave good correspondence 
between estimated allele frequencies in pooled samples and genotypes verified by 
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analysis of the single samples. A similar approach applied to a different application 
is to analyze the Opsin color vision gene array. Within tins gene there are several 
copies on the same chromosome. The tandem array is believed to consist of a single 
red (L) gene, and one to seven copies of the green (M) gene, as a result of uneven 
5 intergenic crossing-over. The red and green genes differ by fifteen amino acid 

residues, and 7 of these 15 are responsible for the absorption spectra. Quantification 
of the gene ratio is demonstrated in FIG. 36. By measuring the area under the peaks 
(alleles), quantitative information about DNA copies entering the PGR reaction is 
obtained. 

10 The same concept has also been applied for the analysis of expression 

differences at mRNA level. By use of heterozygous intra exonic SNPs differences 
in liiRNA expression can be observed and would give similar electropherogram as in 
FIG. 36. 

Allelic Imbalance 

1 5 Study of human cancers has indicated a variety of inactivation of multiple 

tumor suppressor genes for tumor development and/or progression cancer. Loss of 
genes or its function can occur through mutation of one allele together with the loss 
of the second allele, or through homozygous deletion of the gene. To determine 
alleUc imbalance (alleUc loss/reduction or allelic gain) in tumor samples, use of 

20 polymorphic microsatellite markers at gene loci of interest has traditionally been 
used. A disadvantage of using microsatellites is that they are rare (compared to 
SNPs) and generally located in non-coding regions of the genome. To circumvent 
the use of microsatellites we have appUed regular SNP analysis by DCE as a means 
for detemiination allelic imbalance. 



25 



While this invention has hc&a particularly shown and described with 
references to preferred embodiments thereof, it will be understood by those skilled 



wo 2005/047541 



PCT/IB2003/005304 



-78- 

in the art that various changes in form and details may be made therein without 
departing from tiie scope of the invention encompassed by the appended claims. 
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CLAIMS 

n 

What is claimed is: 



1 . A method for separating nucleic acids comprising electrophoresing a sample 
applied to a gel electrophoresis matrix in a capillary^ wherein during 

5 electrophoresis, the temperature of the matrix is cycled at least two times 

between a high and low temperature. 

2. The method of Claim 1, wherein the nucleic acids to be separated are DNA 
fragments comprising one or more polymorphic sites. 

3 . The method of Claim 2, wherein allelic variants at the one or more 
10 polymorphic sites are separated. 

4. The method of Claim 1, wherein the temperature is initially at a high 
temperature and the first cycle is from a higli temperature to a low 
temperature. 

5. The method of Claim 1, wherein the hi^ temperature and/or low 
15 temperature is different during successive cycles. 

6. The method of Claim 1, wherein the temperature is cycled from about 2 to 60 
times. 



7. 



The method of Claim 4, wherein the temperature is cycled about 20 times. 
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8. The metliod of Claim 1 , wherein the high temperature is about 3 ""C higher 
than the low temperature. 

9. The method of Claim 1 , wherein the temperature is between about 2 °C and 
about 15 ""C higher than the lower temperature. 

5 10. The method of Claim 1 , wherein the higher temperature is between about 3 
""C and about ICC higher than the lower temperature. 

1 1 . The method of Claim 1 , wherein the high temperature is less than about 80 
^C. 

12. The method of Claim 1, wherein the low temperature is about 40 °C. 

10 13. The method of Claim 1 , wherein the high temperature is between 50 °C and 
75 °C. 

14. The method of Claim 1, wherein the low temperatm*e is betv^^een 40 °C and 
50 ^G. 

1 5 . The method of Claim 1 , further comprising detecting dsDNA after 
15 electrophoresis. 

1 6. The method of Claim 13, wherein, after the desired number of temperature 
cycles have been completed, the temperature of the gel matrix is such that 
DNA remains double-stranded. 
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17. The method of Claim 1, wherein the temperature oscillations are ramped to 
provide optimal separation of the alleles. 

18. A method for estimating allele frequency comprising: 

electrophoresing a sample applied to a capillary gel 
5 electrophoresis matrix, wherein during electrophoresis, the 

temperature of the matrix is cycled at least two times, wherein one 
cycle is from a high temperature to a low temperature or from a low 
temperature to a high temperature, thereby separating DNA 
molecules in the sample; and 

1 0 quantifying the variant sequences of the separated DNA 

molecules 

thereby providing an estunate of tlie allele frequency for each variant DNA 
molecule. 

19. The method of Claim 18, ftirther comprising detecting dsDNA after 
15 electrophoresis. 

20. The method of Claim 19, wherein, after the desired number of temperature 
cycles have been completed, the temperature of the g:el matrix is such that 
DNA remains double-stranded. 

21 . A method for detecting a microhaplotype comprising separating DNA 

20 fragments comprising a sequence comprising two or more polymorphic sites 

of tlie microhaplotype, wherein the fragments are separated by capillary 
electrophoresis perfomied with two or more temperature oscillations 
between a high and a low temperature. 
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The method of Claim 21, wherein the temperature oscillations are ramped to 
provide optimal separation of the microhaplotype. 
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FIG. 4 
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